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Abstract 


We study lepton mixing patterns which can be derived from the A 5 family symme¬ 
try and generalized CP. We find five phenomenologically interesting mixing patterns 
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for which one column of the PMNS matrix is 


5+V5 1 1 ^1 

’ Vs+Cs’ Vs+Cs 


(the first col¬ 


umn of the golden ratio mixing), (y ^ .iG , —J==, , ^ Y' (the second column of the 

yb—vb V5—v5 

golden ratio mixing), (1,1, /^/^ or (\/5-|-l, —2, \/5 —1)'^/4. The three lepton mixing 

angles are determined in terms of a single real parameter 0 , and agreement with exper¬ 
imental data can be achieved for certain values of 6 . The Dirac CP violating phase is 
predicted to be trivial or maximal while Majorana phases are trivial. We construct a 
supersymmetric model based on A 5 family symmetry and generalized CP. The lepton 
mixing is exactly the golden ratio pattern at leading order, and the mixing patterns of 
case III and case IV are reproduced after higher order corrections are considered. 
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1 Introduction 


In the standard three flavor neutrino oscillation paradigm, lepton flavor mixing is de¬ 
scribed by the so-called Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix Upmns which 
is a 3 X 3 unitary matrix [^. Upmns contains three mixing angles 6^12, 6^13, 623 and one Dirac 
CP violating phase 6cp- There are two more Majorana CP phases if neutrinos are Majorana 
particles. With the measurement of the last mixing angle 633 by Daya Bay |^, RENO and 
Double Chooz , all three lepton mixing angles have been measured with good accuracy in 
neutrino oscillation experiments [s]^. Recently T2K has reported a slight preference for 6cp 
close to 37r/2 [^, when the data are combined with the measurements of the reactor experi¬ 
ments. The present global fit to neutrino data also indicates nontrivial values of Scp [5]-[^. 
However, the values of the both Majorana phases are unknown so far. Search for leptonic 
CP violation via the determination of Scp is one of the major goals of future long-baseline 
experiments such as the proposed LBNE [^, LBNO 


10 and HyperKamiokande 11 


Now it is established that both neutrino and charged lepton mass matrices have residual 
flavor symmetries determined by lepton flavor mixing, and vice versa residual flavor symme¬ 
tries in the mass matrices can determine the lepton mixing matrix up to Majorana phases 
and permutations of rows and columns [^. Inspired by the fact, it is assumed that the 
residual flavor symmetries arise from a underlying flavor symmetry group Gf which is usu¬ 
ally chosen to be a finite and non-abelian subgroup of U{3). In the past years, much effort 
has been devoted to the discussion of lepton flavor mixing from a discrete flavor symmetry 
Gf and its breaking 13 . It is surprising that the mixing patterns achievable in this way are 


quite limited, the PMNS matrix can only be of the trimaximal form to accommodate the 


experimental data and the Dirac phase is trivial 14 


Beside residual flavor symmetries, neutrino and charged lepton mass matrices have resid- 

Analogous to residual flavor symmetries, residual CP symmetries 


ual CP symmetries 15,16 


also impose strong constraints on the mass matrices and therefore allow us to reconstruct 
the lepton mixing matrix 15 . A simple example is the well-known /i — r reflection sym¬ 
metry 17 which predicts maximal atmospheric mixing angle 623 and maximal Dirac CP 


phase. It is natural to conjecture that there is a CP symmetry Hcp (also called generalized 
CP symmetry) at high energy scale, which is broken down to the residual CP symmetries 
at low energy. Note that the effects of CP symmetry on the fermion mass matrix have been 


discussed several decades ago 18,19 


Recently it is proposed to predict the lepton mixing angles and CP phases by combining 
a discrete flavor symmetry Gf with a CP symmetry Hcp 20,21 . Hcp has to be compatible 


with Gf such that the possible forms of the CP transformations are strongly constrained. 
It has been proved that the mathematical structure of the group comprising Gf and Hcp 
is in general a semi-direct product G/ x Hcp 


20 


In this framework, the flavor symmetry 
Gf is broken down to different abelian subgroups G^ and Gi in the neutrino and charged 
lepton sectors respectively, and Hqp is broken into residual CP symmetry Hqp and 
respectively. The mismatch between the remnant symmetries Gi, x H^p and Gi x iP^p 
generates the PMNS matrix. Neutrinos are generically assumed to be Majorana particles. 
As a consequence, Gy can only be a or Z2 subgroup of Gf. In the case that Gy = K4 
and Gi is capable of distinguishing the three generations (i.e.,G/ can not be smaller than 
Z3), all lepton mixing parameters including the Majorana phases would be completely fixed 
by residual symmetries once the CP symmetry is considered. In this way, both Dirac and 
Majorana CP violating phases are found to be conserved in the context of A(6n^) family 
symmetry combined with generalized CP [^. Recently a bottom up analysis of the remnant 
K4 flavor symmetry and CP symmetry in the neutrino sector has been performed 16 . On 
the other hand, if Gy = Z2 and a CP symmetry is preserved in the neutrino sector, only 
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one column of the PMNS matrix can be fixed and all lepton mixing parameters depend on 
one single real parameter 9. Along this line, the family symmetries A 4 [^, 5*4 20,24 -28 , 
T' 1^, A(48) 1^, A(96) 1^, A(3n^) and A( 6 ? 7 ,^) 1^,33 which are combined with the 


corresponding generalized CP symmetries have been investigated already. It is found that CP 
phases can only be trivial or maximal in simple family symmetries A 4 23 and 5*4 20,24-28 
while A(48) and A(96) (also A(3n^) and A( 6 n^) [^[^) family symmetries admit 


mixing patterns in which all CP phases nontrivially depend on the parameter 9. In addition, 
some models with both flavor and CP symmetries have been constructed 24-28,30 . Last 


but not least, if remnant symmetries in the neutrino and charged lepton sectors are 7^4 x H^p 
and Z 2 X respectively, then the PMNS matrix is also predicted in terms of the parameter 
9 and one row instead of one column would be hxed 


28,33 


It is known that the flavor symmetry group should be of the von Dyck type 
hnite von Dyck groups include S 3 , A 4 , 54 , A 5 and dihedral groups 


34 . The 


35 . Since S 3 and 


dihedral groups don’t have irreducible three dimensional representations, they are not suit¬ 
able as flavor symmetry otherwise two mixing angles would vanish. The phenomenological 
consequences of A 4 and 5*4 flavor symmetries combined with generalized CP have been stud¬ 
ied 


20,23 


In the present work, we shall investigate the A 3 flavor symmetry and CP 
symmetry. We shall perform a model independent analysis of possible lepton flavor mixing 
obtained from breaking of the original symmetry A3 x Hqp- We hnd hve phenomenologi¬ 
cally interesting mixing patterns summarized in Table The three mixing angles turn out 
to depend on only one free parameter 9 and good agreement with their measured values can 
be achieved for certain values of 9, the Dirac CP phase is conserved or maximal and the 
Majorana CP phases are trivial. Furthermore, we construct a model based on A 3 x Hcp- 
The lepton mixing is exactly the golden ratio (GR) texture at leading order (LO). A non-zero 
9 i 3 is generated by the next-to-leading-order (NLO) corrections, and the mixing patterns of 
cases III and IV discussed in the model independent analysis are generated. 

The layout of the rest of this paper is as follows. In sectionthe physical CP transforma¬ 
tions compatible with the A 3 family symmetry are found. In section we perform a model 
independent analysis of possible lepton mixing patterns achievable from the underlying sym¬ 
metry group A3 X Hcp- In section we present our A3 x Hcp model, the LO structure, 
vacuum alignment and the NLO corrections are discussed. Section concludes the paper. 
In Appendix]^ we review the group theory of A 3 and the Clebsch-Gordan coefficients in our 
working basis are reported. In Appendix we present the possible mixing patterns arising 
from the A 3 flavor symmetry without CP symmetry, where the residual flavor symmetry in 
the neutrino sector is either Klein or Z 2 subgroup of A 3 . Compared with section]^ we see 
that generalized CP is really a powerful method of predicting CP phases as well as lepton 
mixing angles. 


2 Approach 

Both family symmetry and CP symmetry acts on the flavor space in a non-trivial way, and 
the interplay between them should be carefully treated. In order to consistently combine a 
family symmetry Gf with a CP symmetry which is represented by unitary CP transformation 
matrix X, X must be related to an automorphism u : Gf ^ Gf. To be precise, the CP 


transformation X should be a solution to the consistency equation 20,21 

A>-(g)A'-‘=p(u(g)), VgeG/, 


( 2 , 1 ) 


where p is a representation of Gf with p : G ^ GL{N,C), and it is generally reducible. 
We can easily check that the automorphism associated with p{h)X for any h G G/ is an 
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composition of u and an inner automorphism fih ■ 9 ^ hgh ^ with h,g E Gj 28,31 . 


Therefore the effects of inner automorphism can be easily included, and it is equivalent to 
a family symmetry transformation. As a consequence, we could hrstly focus on the the 
outer automorphism of Gf. Furthermore, it has been shown that u has to a class-inverting 
automorphism for A to be a physical CP transformation 


36 


In other words, u should 
map each irreducible representation r of Gf into its own complex conjugate. Hence the 


consistency condition in Eq. ( 2.1[ ) takes a more restricted form: 

XrPrig)^;:' = Pr{u(g)) , yg€Gj, 


( 2 , 2 ) 


where the subscript “r” refers to the representation space acted on. The CP transformation 


X in Eq. (2.1) is given by the direct sum of the Xr corresponding to the particle content of 


the model. Notice that the consistency conditions of Eq. (2.2) can also be derived from the 


requirement that the Lagrangian is invariant under both CP symmetry and flavor symme¬ 


try 37 


In the present work, we are interested in the family symmetry Gf = A^. The group 
theory of A^, its representation and all the Clebsch-Gordan coefficients are reported in 
Appendix The structure of the automorphism group of A^ is quite simple and is very 
clear in mathematica. 


2(^5) = Zi, Aut(A5) = Fs, 

Inn(A 5 ) = A 5 , Out(A 5 ) = Z 2 , (2-3) 

where Z{A^), Aut(A 5 ), Inn(A 5 ) and Out(A 5 ) denote the center, automorphism group, inner 
automorphism group and outer automorphism group of A 5 respectively. We see that the 
outer automorphism group of A 5 is isomorphic to Z 2 . Consequently there is only one non¬ 
trivial outer automorphism u with 

T^{ST^Y. (2.4) 

The order of u is really 2, i.e., = id, where id represents the trivial automorphism 

id{g) = g,\/g E A 5 . One can straightforwardly check that u acts on the A 5 conjugacy classes 
as follows 

ICi ^ IGi, I 5 G 2 ^ ISC's, 2 OC 3 ^ 2 OC 3 , I 2 C 5 M I 2 C 5 . (2.5) 

It interchanges the classes I 2 C 5 and I 2 C 5 . Since the inverse of each A 5 conjugacy class is 
equal to itself, u is not a class-inverting automorphism, and the corresponding CP trans¬ 
formation is unphysical. In terms of representations, the two different three-dimensional 
irreducible representations 3 and 3' are exchanged not mapped into their conjugate under 
the action of u. The generalized CP symmetry related with u can only be consistently de- 
hned if helds transforming as 3 and 3' are absent in a model. As a result, we conclude that 
only the CP transformation associated with the trivial outer automorphism (i.e., the inner 
automorphism) can be compatibly imposed on the theory with A 5 family symmetry. 

Now we consider the representative inner automorphism {.S,T) —)■ {S,T^). 

The corresponding generalized CP transformation X^ is fixed by the consistency equations: 

= pjs), 

X,“p;(T)(Af)-‘=p,(T‘‘). (2.6) 

From the representation matrices given in Appendix we see that for any representation 

p:(S)=MS), p:(T)=MT^). (2.7) 
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( 2 . 8 ) 


Therefore is an identity matrix np to an overall phase, i.e., 

= 1. 

Inclnding the contribntion of the remaining inner antomorphisms in the manner stated below 


Eq. (2.1), the most general CP transformation consistent with family symmetry is of the 
form 

= Prig)X^ = Pr{g), g G . (2.9) 

This means that the generalized CP transformation consistent with A^ is of the same form as 
the family gronp transformation in onr working basis while they act on the a held mnltiplet 
in different ways: (p{x) \— > pr{g)(p{x), g E A^ versns (p{x) \— > Xj.ip*{xp) = pj.{g)(p*{xp), 
where xp = {t, —x). 

In this work, the phenomenological implications of A^ family symmetry combined with 
the generalized CP symmetry wonld be investigated in a systematical and comprehensive 
way. The parent symmetry is A^ x Hcp at high energy scale, where the element of Hcp is 


the CP transformation compatible with A^ and its explicit form is given by Eq. (2.9). In 


this setnp, lepton mixing can be predicted from A^ x Hqp breaking into different remnant 
symmetries Gi x H]jp and Gy x H^p in the charged lepton and nentrino masses respectively, 
where G;, Gu and H\jp, H^p denote residnal family symmetries and residnal CP symmetries 
respectively. It is notable that the predictions for the lepton havor mixing only depend on 
the assnmed symmetry breaking patterns and are independent of the details of a specihc 
implementation scheme, snch as the possible additional symmetries of the model and the 
involved flavon helds and their assignments etc. In practice, the three generations of left- 
handed leptons donblets are embedded into the faithfnl three-dimensional representation 
3 of A 5 . Since 3' is related to 3 by the onter antomorphism u, the resnlts wonld be the 
same and no additional resnlts would be found, if we assign the three left-handed leptons to 
the representation 3' instead. The requirement that Gi x H^^p is preserved by the charged 
lepton mass term implies that the hermitian combination m\mi must be invariant under the 
remnant symmetry Gi x H^p, i.e.. 


pl{gi)m]mips{gi) = m\mi, 

Xj^mlmiXia = 


gi £ Gi 
Xis e m 


CP ) 


(2.10a) 

(2.10b) 


where the mass matrix mi is dehned in the convention Ipmili. Once Gi and H^ p are specihed , 


the most general form of mlmi can be straightforwardly constructed from Eqs. (2.10a 


2.10b). 


In the present work, we shall assume neutrinos are Majorana particles. In the same fashion, 
requiring that Gi, x H^p is a symmetry of the neutrino mass matrix m,^ implies that 
should be invariant under the action of G,, x H, 


CP-: 


pl{gu)muP3igt,) = m^, g^EGy, 

Xl^m^X^s = ml, X^s E HI 


CP y 


(2.11a) 

(2.11b) 


which allow us to derive the explicit form of m„. Since both remnant family symmetry 
and remnant CP symmetries are still preserved after symmetry breaking, they should be 


compatible with each other. That is to say consistency equation similar to Eq. (2.2) has to 
be fulhlled. 


Xup*{gu,)xp = p{g^^), 

Xip*{giDX^^ = p(gi^), 


gvitg^j ^ Gi/, 
giiigu ^ Gi. 


(2.12a) 

(2.12b) 


The prediction for the PMNS matrix can be obtained by further diagonalizing the recon- 

i\mi and Please see Ref. [1^ for an alternative way of directly 


structed mass matrices m]mi and m,. 
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extracting the PMNS matrix from the representation matrices of the remnant symmetries 
without resorting to the mass matrices. As the order of neutrino and charged lepton masses 
is indeterminate in our framework, it is only possible to determine the PMNS matrix up to 
independent row and column permutations. 


From the remnant symmetry invariant conditions of Eqs. (2.10a, 2.10b), we can see that 
Xiy and pr{gi)Xir with gi G Gi lead to the same constraint on mjm^. ^Furthermore, the 
residual CP transformation X/r should be a symmetric matrix otherwise the charged lepton 
masses would be restricted to be partially degenerate 15,28 . The same comments apply to 
Xijj- and pj-i^gif^Xif 


with g^ G Gy. Notice that the same result for PMNS matrix would be 
obtained (M,28,^, if a pair of subgroups is conjugated to the pair of subgroups 


{Gi,Gy} under an element of A^, i.e.. 


G[ = gGig 


-1 


G'y — gGyg 


-1 


g E A 5 


(2.13) 


The reason is that remnant CP symmetries determined by restricted consistency condition 


of Eqs. (2.12a, 2.12b) are strongly correlated in the two cases such that lepton mass matrices 
{m'lm'i, m'y} for the new primed residual symmetry are related to rriy} by a similarity 

In this way, it is sufficient to only discuss the independent 


transformation psig) [23,28,31 


pairs of {Gi,Gy} which are not related by group conjugation and subsequently all possible 
residual CP compatible with the residual family symmetry should be included . 


3 Lepton mixing from remnant symmetries of ^5 Hcp 

Neutrino are assumed to be Majorana particles here, therefore the remnant flavor sym¬ 
metry Gy must be a Klein four 7^4 = Z 2 x Z 2 subgroup or a single Z 2 subgroup of A 5 . Gi can 
be any abelian subgroups of A^ with order equal or greater than 3. A complete or partial 
degeneracy of the charged lepton mass spectrum would be produced if Gi had a non-abelian 
character. In the case of Gy = the lepton mixing matrix Upmns is fully determined by 
the mismatch between the remnant family symmetry Gi and Gy. As shown in Appendix [B| 
UpMNS can take four possible forms such as the golden ratio mixing, democratic mixing 
and so on. However, none of them is compatible with experimental data. Then we turn to 
the scenario of Gy = Z 2 . With this setting, Upmns is partially constrained, and only one 
column of the lepton mixing matrix is hxed up to reordering and rephasing of the elements. 
The explicit forms of the hxed column vectors for all the independent residual havor sym¬ 
metries are summarized in TableWe hnd that four cases are viable: (G;, Gy) = (Zf’, 2’f), 
(Zj’, and \ Zf) lead to the mixing column 

vectors (-^/^, _A=)r ( ( 1 1 1 )T ^nd 


respectively, where k = (1 -|- -\/5)/2 is the golden ratio. The phenomenological implications 
of each case are explored in Appendix [^ and the lepton mixing matrix Upmns turns out 
to depend on two free parameters up to indeterminant Majorana phases. We see that the 
measured values of the three mixing angles can be accommodated very well, but the allowed 
values of Dirac CP phase 5cp scatter in a quite large range. Furthermore, the breaking 
patterns with {GpGy) = (Z 2 ,iC 4 ) are studied as well. Accordingly a row of the lepton 
mixing matrix Upmns is determined to be f (k, — 1) or (1,0,0) which are not in the 

experimentally preferred regions. 

In order to be able to predict the values of CP phases, we extend the A 5 family symmetry 
to include the generalized CP. In the following, we shall perform a thorough analysis of lepton 
mixing patterns for the possible residual symmetries Gi x H^p and Gy x H^p in the charged 
lepton and neutrino sectors, where the remnant family symmetries Gi and Gy would be 
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CP 


restricted to the four viable cases listed in Table and the remnant CP symmetries H[ 
and H^p are determined by consistency condition of Eqs. ( 2 . 12 a[ 2 . 12 b ). In this setup, Upmns 
as well as all mixing angles and all CP phases generically depend on a free parameters 6 
whose value can be fixed by the measured value of 613 . As a consequence, all observables 
are strongly correlated. For the concerned A 5 family symmetry, the Dirac phase would be 
predicted to be trivial or maximal while both Majorana phases are trivial after generalized 
CP symmetry is imposed. In order to evaluate how well the predicted mixing patterns agree 
with the experimental data on mixing angles, we shall perform a usual analysis which 
uses the global fit results of Ref. [^. We begin to discuss all possible cases one by one. 


3.1 G] — , Gjy 


= zi 


i — ^5 ■> 

In this case, the parent symmetry A 5 xi Hcp is broken down to xi H^p and Z 2 x 
Hqp subgroups in the charged lepton and neutrino sectors, respectively. The residual CP 
symmetry H\jp must be consistent with the residual flavor symmetry Z'^ in the charged 
lepton sector. That is to say the element Xir of Hjjp should fulfill the consistency equation 


of Eq. (2.12b), 


WrP:(T)X-i=Pr(^'), 


9 ' e Zj 


(3.1) 


Then we hnd only 10 choices out of the 60 CP transformations of Hcp listed in Eq. (2.9) 
are acceptable 


{p,(i),pAT),pAP)<p.(P).pAP),p,(sr-sPS),pMT'‘s)Vs), 
Pr(PSPSPS),Pr('PSp-SPS),Pr(SPSPS )} . 


(3.2) 


As shown in Eq. (2.10a), the residnal family symmetry impose the following constraint 
on the charged lepton mass matrix: 


P3(T)mJmiP3(T) ^ m]mt . 


(3.3) 


In our working basis, the representation matrix of the generator T is diagonal with p^{T) = 
diag(l, (Us, (Ug). Consequently the hermitian combination m\mi of charged lepton mass matrix 


is also diagonal, i.e.. 


m\mi = diag (rUg, m^) , 


(3.4) 


where rUe, and ruy represent the electron, muon and tau masses respectively. Fur¬ 


thermore, we can check that the remnant CP invariant condition of Eq. (2.10b) is auto¬ 


matically satisfied for Xi^ = priV)) PriT), priT"^)) PriT^)) PriT^)■ However, the mass de¬ 
generacy = rriT- arises for the remaining values Xiy. = p^ST'^ST^S), pr{(T^SYT^S), 
PriT^ST^ST^S), priT^ST^ST^S), pr{ST^ST^S). The reason is that all remnant CP trans¬ 
formations except pr{T^ST"^ST^S) are not symmetric. Generally speaking, any remnant CP 
transformation must be a symmetric matrix to avoid degenerate masses 15,28 . This case 
is obviously not viable, and will be disregarded hereafter. 

Now we turn to the neutrino sector. The residual CP transformations X, 
specified by the consistency condition: 


of H^p is 


A'„p;(S).Y-‘ = p,(S ), 


(3.5) 


which can be easily obtained by applying the general consistency condition of Eq. (2.12a). 
We see that the CP transformation X, 


commutes with flavor symmetry transformation 
Pr{S), and therefore remnant symmetry is Zf x H^p in the neutrino sector in this case. 
Notice that the semi-direct product structure between residual flavor and CP symmetries 
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generally reduces to a direct product if the residual flavor symmetry is a Z 2 subgroup 23,24 
It is easy to check that X^r can only take 4 possible values, 


H'^p = {Pr{l),Pr{S),p,{T^ST^ST^),p,{T^ST^ST^S)}. (3.6) 

The neutrino mass matrix rriv respects the residual symmetry Zf x satisfying 

pl{S)mup^{S) = rriu , 

Xjgm.X.a = <, X.3 e H^cp • (3-7) 


We hud that the most general neutrino mass matrix invariant under the residual family 
symmetry Zf, takes the following form 

/100\ . /-2V2 33\ /2 0 0\ /O -^/2 V2\ 

m^ = a 001+^ 3 0^2 +7 0 3 -1 +5 -^2 -2k 0 , 

\010/ ^2\3 V 2 0 J \0 -1 3 J \v^ 0 2k J 

(3.8) 

where a, f3, 7 and 6 are generally complex parameters, and they are further constrained to 
be real or pure imaginary by residual CP. This neutrino mass matrix rrii, can be simplihed 
into a quite simple form by performing a golden ratio transformation. 


= Uaj^m^GR = 


a — {3 k — 1)(3 + 27 0 

0 a + {3 k — 2)/3 + 27 

0 2 v/ 2 T^ 5 


0 


2y/2^ 5 
/? + 47 


-a 


where 


Ugr — 


V Vs 


2+5 f 


\ V 2^/5A 


^ 0 \ 


1 
Vs 

K. 

2+5 

K 

2+5 


1 

V2 

I 2 / 


(3,9) 


(3.10) 


is the golden ratio mixing pattern 38 which can be naturally derived in models pl9] . The 


neutrino mass matrix m'^ is further diagonalized by a unitary rotation f/(, in the (2,3)-plane, 

U'Jm'JJl = diag(mi, m 2 , m 3 ). (3.11) 


The next step is to explore the constraint of remnant CP on rriy. Two different phenomeno¬ 
logical predictions arise for the four possibe X^r shown in Eq. (3.6), as pr{S)Xyr and Xyr 
lead to the same predictions. 


(I) X,r = Pr(l),Pr(5) 

Obviously we have m,^ = m* such that all the four parameters a, /3, 7 and 5 are real. 
As a consequence, the neutrino mass matrix m'^ is a real symmetric matrix. The unitary 
transformation f/' is of the form: 


/I 0 0 \ 

0 COS0 sin0 A:^. (3.12) 

y 0 — sin 9 cos 9 J 

where K^, is a diagonal phase matrix with elements equal to ±1 or ±i which makes the 
neutrino masses mi_ 2,3 positive. The effect of K^, is a possible change of the Majorana 
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phases by vr, and it would be omitted hereinafter for the other cases. The parameter 6 is 
given by 

4:\/2 + K 6 , , 

^ (3.13) 


tan 20 = 


2{a — 7 ) + (3k — l)j3 


The light neutrino mass eigenvalues are 

mi = \a — (3k — l)/3 + 27 I , 


mo = - 
2 


m 3 = 


3(k - l)(3 + 67 + 
3(k — l)/d + 67 — 


2(a — 7 ) + (3k — l)/d 


cos 29 

2{a — 7 ) + (3k — l)/3 


cos 20 


(3.14) 


Given the diagonal charged lepton mass matrix, the lepton mixing matrix takes the form 




UpMNS — UgrUI — 




2 VbK 




V 


K. 


(3.15) 


K O sin 6 / K I COS 9 

^COS 0 -^ + ^ j 


2V5k Y 2 V 5 \/2 Y ' 

One can straightforwardly extract the lepton mixing angles and CP phases as follows, 

•2/) 3-K.2. . 2n 1 + COS 20 

sm 013 = -sm 0 , sm 0 i 2 = 


sin^ 003 = V 


Y2 + K sin 20 


3k — 1 + (k — 1) cos 20 ’ 


3 + 2k + cos 20 ’ 
sin(5cp = sin 0)21 = sin 031 = 0 , (3.16) 


where 6 cp is the Dirac CP phase, a 2 i and 031 are the Majorana CP phases in the standard 
parameterization [^. There is no CP violation in this case as the neutrino mass matrix 
is real. Expressing 0 in terms of 0 i 3 , correlations among the three mixing angles follow 
immediately. 


sim 010 = 


K 


2 + K 


tan^ 013, 


sin^ 023 = - ± K tan 0 i 3 •^l — (1 + k) tan^ 0 i 3 . 


(3.17) 


For the measured reactor mixing angles sin^ 0 i 3 ~ 0.0234 [^, we have sin^ 023 — 0.258 or 
0.742 which is outside of the experimentally favored 3 ct region although sin^ 0 i 2 — 0.259 
is acceptable. As a consequence, this mixing pattern isn’t viable. This point remains even 
after permutation of rows and columns is considered. 


(II) = Pr{T^ST^ST^),pr{T^ST^ST^S) 


Solving the residual CP invariant condition in Eq. (3.7), we hnd a, 13 and 7 are real while 6 
is pure imaginary. The unitary diagonalization matrix U'^ is 

u; = 


1 0 0 

0 cos 0 sin 0 
0 — i sin 0 i cos 0 


(3.18) 


where the diagonal matrix K^, multiplied from the right-hand side has been omitted, and 
the rotation angle 0 fulhlls 

4i\/2 + K 6 
3(k — 1)(3 + 67 


tan 20 


(3.19) 
































The three neutrino masses are given by 


mi = \a — (3k — !)/? + 27 I , 
m 2 


m 3 = 


2a + (3k — 1)13 — 27 + 
2a + (3k - 1)13 - 27 - 


3((k-1)/? + 27) 


cos 29 

3 ((k - l)/3 + 27 ) 


cos 29 


(3.20) 


All the four parameters a, (3, 7 and 6 are involved in the three neutrino masses. As a result, 
the measured mass squared differences 6 m‘^ = m\ — and Am^ = m^ — (m^ + m |)/2 
can be easily accommodated [^, the absolute neutrino mass scale can not be hxed, and the 
neutrino mass spectrum can be either normal ordering (NO) or inverted ordering (10). The 
PMNS matrix takes the following form: 


U. 


PMNS 


V Vs 


2^5 f 


\ V 


COS 9 

V5«: 


^COS 0 +^ 


sin 9 \ 

VOK 


O icosO 

2V5 V 2 




(3.21) 




Note that the hrst column vector of this mixing pattern coincides with the hrst column of 
the GR mixing. The lepton mixing angles and CP phases can be read out a^ 


3 — K 


sin^ 013 = 


sin^ 023 = X , 


sin^ 0 , 


sin 012 = 


1 + cos 20 


3 + 2k + cos 20 ’ 
|sin Scp\ = 1 , sin 0:21 = sin a^i = 0 . 


(3.22) 


Here we present the absolute value of sirnfcp, since the sign of simf^p depends on the 
ordering of rows and columns. We see that both atmospheric angle 023 and Dirac CP 
phase 6 cp are maximal while Majorana phases are conserved. Given the weak evidence of 
5cp ~ 37 r /2 from T2K [^, this pattern is slightly preferred. The prediction of maximal 
Dirac CP can be tested by next generation long-baseline neutrino oscillation experiments 
such as the proposed LBNE [^, LBNO 


10 and HyperKamiokande 11 , which aim to search 


for leptonic CP violation. Moreover, the correlation between 0 i 3 and 0 i 2 is of the same 
form as that of case I, and it is plotted in Fig. [T] The results of the analysis are 
reported in Table We see that the experimental data on lepton mixing angles can be 
accommodated very well. Notice that the solar mixing angle 0 i 2 is predicted to be around the 
present 3a lower bound. As far as we known, the JUNO experiment can measure 0 i 2 with 
high accuracy 44 . If signihcant deviations sin^ 0 i 2 from 0.259 was detected, this mixing 
pattern would be excluded. It is well-known that leptonic CP phases can play a crucial 
role in the rare process neutrinoless double beta {{(3(3)oi,—) decay. The dependence of the 
(/d/3)ojy—decay amplitude on the neutrino mixing parameters is characterized by the effective 
Majorana mass |mee| with the dehnition: 

\mee\ = \mi COS^ 012 COS^ 013 J- m2 siu^ 012 COS^ 0136 *"^^ J- m 3 siu^ 5)^3e*("3i-25cp) I _ (3.23) 


For the predicted mixing parameters in Fq. (3.22), we have 


|mee| = Umi-|-K ^A:2^2r2 cos^ 0-|-K ^fc3m3sin^0| , 
V5 


(3.24) 


^In the case of sin 20 = 0, either 0 i 2 or 0 i 3 vanishes, consequently the value of Scp can not be determined 
uniquely. 
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Figure 1: The correlation between sin^ 612 and sin 6*13 (left panel) and the allowed values of the effective mass 
\mee\ (right panel) in case II. On the left panel, the best fitting value is labelled with a red pentagram, and 
the points for 0 = 0, tt/G, tt/S and 7 r /2 are marked with a cross to guide the eye. The Icr and Scr ranges of the 
mixing angles are taken from Ref. . On the right panel, the orange and green bands denote the 3 ct regions 
for normal ordering and inverted ordering mass spectrum respectively. The red and purple areas are the 
predictions for the lepton mixing matrix in Eq. (3.21). The present most strict bound \mee\ < (0.120 — 0.250) 
eV from EXO-200 40 41 combined with KamLAND-ZEN is represented by the horizontal dashed line, 
and the upper limit on rrimin from the latest Planck result mi + m 2 + m 3 < 0.230 eV at 95% confidence 
level 


43 


from the latest Planck result mi 
is shown by vertical dashed line. 


where ^ 2,^3 = ±1 originates from the ambiguity of the matrix K^. The prediction for the 
effective mass |mee| with respect to the lightest neutrino mass is shown Fig. We hnd that 
\mee\ is close to 0.022eV or the upper bound 0.045eV in case of 10 neutrino mass spectrum, 
which are within the future sensitivity of planned (/d/3)oi/—decay experiments. However, in 
case of NO spectrum, |mee| strongly depends on lightest neutrino mass rrimin, and it can 
even be approximately vanishing for particular value of rrimin- 


3.2 


Gi = Zi ,G, = Z, 


T^ST'^ST^ 


The charged lepton sector preserves the same remnant symmetry Zj x 


H\jp as that 


discussed in section [TT| Therefore the charged lepton mass is subject to the same constraint, 
and m\mi should be diagonal as well. In neutrino sector, the residual CP symmetry H'^p 
has to be compatible with the residual family symmetry Gi, = , i.e., 


= pZt-^st^st 


(3.25) 


It is easy to check that only 4 generalized CP transformations are acceptable, 

H^cp = {pril),Pr{S),pZT^ST^ST^),pZT^ST^ST^S)} . (3.26) 

Straightforward calculations demonstrate that the most general neutrino mass matrix in¬ 
variant under form 

/10 0\ g /- 2 V 2 33\ /2 0 0\ / O V2k -V2k\ 

m^ = a 001+^ 3 0^2 +7 0 3 -1 +5 y/2K -2 0 h 

\0 10 / 3 ^ 20 / \0-13/ \-y/2K 0 2 / 

(3.27) 


10 































Analytic expression 

Best fitting 


sin^ 6 »i3 

sin^ 6*12 

sin^ 023 

Obf 

'X.min 

sin^ 013 

sin^ 012 

sin^ 023 

II 

10 

^ sin^ 6 

2 cos^ 9 
3+2k+cos 29 

1 

2 

0.295 

8.468 

0.0234 

0.259 

0.5 

NO 

1 

2 

0.292 

11.88 

0.0229 

0.259 

0.5 

III 

10 

^ sin^ 6 

V5 

4:-2k 

5 — 2k+cos 29 

1 \/3 —K sin 29 

0.182 

4.851 

0.0236 

0.283 

0.404 (023 < 45°) 

2.958 

3.165 

0.0240 

0.283 

0.597 (023 > 45°) 

NO 

2 3k—2-\-kcos29 

0.179 

4.087 

0.0230 

0.283 

0.406 (023 < 45°) 

2.965 

24.88 

0.0224 

0.283 

0.593 (023 > 45°) 

IV 

10 

1 

2 

0.183 

2.232 

0.0241 

0.283 

0.5 

NO 

1 

2 

0.181 

5.802 

0.0235 

0.283 

0.5 

V 

10 

1 —sin 2Q 

3 

1 

2+sin 29 

1 

2 

0.976 

3.987 

0.0238 

0.341 

0.5 

NO 

1 

2 

0.973 

7.480 

0.0233 

0.341 

0.5 

VII 

10 

(cos 6 — k sin 0)^ 

(k cos @+sin 9)^ 

(k^ cos 0—sin 

4k. 2 —(cos 9 — k sin 9)^ 

0.286 

1.626 

0.0242 

0.329 

0.486 (023 < 45°) 

K^{cos 9 -\-k sin 9)^ 
4/s^^ — (cos 9 — k sin 0)^ 

0.286 

1.751 

0.0242 

0.329 

0.513 (023 > 45°) 

NO 

4k2 

4k, 2 — (cos 9 — k sin B)‘^ 

{k^ cos 6—sin 9)^ 

4:K^ — {cos 9 — k sin 0)^ 

0.293 

3.503 

0.0229 

0.330 

0.480 (023 < 45°) 

/si^(cos 9-\-k sin 9)'^ 
4k.2 — (cos 9 — k sin 9)^ 

0.282 

6.958 

0.0248 

0.329 

0.5 1 0 (023 > 45°) 


Table 1: Summary of the predictions for the lepton mixing angles and their best fitting values for all viable 
cases in the framework of x Hcp- In case VII, the mixing patterns for 6*23 in the first and second octant 
are related through the exchange of the second and third rows of the PMNS matrix. Notice that all the 
three CP phases are independent of 6 in all cases: Dirac phase is trivial or maximal, and both Majorana 
phases are trivial. 


where the parameters a, /3, 7 and 6 are generically complexes, and they are further con¬ 
strained by the remnant CP. After performing a GR transformation, becomes 

/a-{3 k- l)f3 + 27 0 2 ^/2TK S \ 

^'u = UcifTiMcR = ( 0 a + {3 k - 2)(3 + 2'J 0 . (3.28) 

y 2 a/2 -|- k 5 0 —a — (3 + A'y J 


In the following, we proceed to investigate the constraints imposed by the remnant CP 
transformations shown in Eq. (3.26). The four possible can be divided into two classes. 


(Ill) X,r = Pr(l),Pr(T3AT2,ST3) 

In this case, the residual flavor and residual CP transformations are of the same form. As a 
result, the four parameters a, /?, 7 and 6 are all real. The neutrino mass matrix can be 
diagonalized by a unitary transformation 


cos 6 0 sin 6 
U'„=\ 0 10 

— sin 6 0 cos 6 


with 


tan 29 = 


4a/2 + K 6 


The three neutrino masses are 


m, = - 


2(a-7) - {3k-2)13 ’ 

2{a - 7 ) - (3fi: - 2)13 


(3.29) 


(3.30) 


—3k(3 -b 67 -b 


cos26' 
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(3.31) 


7712 = \a + (3k - 2)(3 + 27 I , 

1 „ „ „ 2 (a - 7 ) - (3k - 2 )/? 


m 3 = 


-3k/3 + 67 


cos 26* 


The absolute neutrino mass scale can not be predicted. Then the PMNS matrix reads 




UpMNS — UgrUI — 


7 sm» \ 


COS 6 


sin 0 


I sin 6 / Av 

V2 Y 2x/5 V2 


COS 0 


\ a/^ 

Note that the second column vector is 


sin 6 

V 2 


Sill (7 / ^ sin 9 I 

“7^ 'i / ^ “T 


cos ^ 

V 2 

cos 0 


(3.32) 


2 vT ^ P2 / 


Tfc’ V V ^J coincides with the 

second column of the GR mixing. The lepton mixing parameters are predicted to be 


■ 2 a ^-2/1 -2/1 4 2k 

sm 6^13 = -G= sm 6 , sm 612 = -, 

^75 5 - 2k + cos 20’ 


sin^ 623 = - - 


a/ 3 — Ksin 26 
3k — 2 + K cos 26 


sin 6 c p = sin a 2 \ = sin 031 = 0 


(3.33) 


We see that 623 deviates from maximal mixing and all the three CP violating phases are 
trivial due to a common CP transformation Pr(l) of fhe charged lepton and neutrino sectors. 
The mixing angles 6 \ 2 i ^13 and 623 only depend on the parameter 0, and they fulhll the 
following relations, 

sin^ 019 cos^ 013 = ^ ^ 


13 = ^—, sin^ 023 = - ± (k - l)tan0i3Yl + (k - 2)tan2 0i3 , (3.34) 


which are plotted in Fig. Obviously the mixing angles can be very close to the their 
measured values for certain values of the parameter 0. The global minimum of the 
function is rather small, as shown in Table The predictions for the effective mass |mee| 
are also displayed in Fig. 

(IV) X,r = Pr{S),Pr{T^ST^ST^S) 

Invariance of the neutrino mass matrix rrii, under the action of these residual CP transfor¬ 
mations implies that a, f3, 7 are real while 6 is pure imaginary. The diagonalization matrix 
of is 


i cos 0 0 i sin 0 
U'„=\ 0 1 0 

— sin 0 0 cos 0 


where 


tan 20 = 


4*72 + K 5 

■3(k/3-27) • 


(3.35) 


(3.36) 


The neutrino masses are given by 

1 


mi = 


2a + (3k — 2)13 -|- 27 -|- 


3(k/3 - 27 ) 


7712 = \a + (3k - 2)(3 -1- 27I , 


m 3 = 


-2a + (3k — 2)13 27 


cos 20 

3(k/? - 27 ) 


cos 20 


(3.37) 
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Figure 2: The correlation among sin^ 0 i 2 , sin^ 6*23 and sin 6*13 (the former three panels) and the allowed 
values of the effective mass |mee| (the last panel) in case III. The global minimum of the function is 
labelled with a red pentagram, and the points for 0 = 0, 7 r/ 6 , tt/S, 7 r/ 2 , 2n/S and 57 r /6 are marked with a 
cross to guide the eye. The la and 3cr ranges of the mixing angles are taken from Ref. [^. In the last panel, 
the orange and green bands denote the 3cr regions for normal ordering and inverted ordering mass spectrum 
respectively. The red and purple areas are the predictions for the lepton mixing matrix in Eq. (3.321. The 
present most strict bound 


rUe 


< 


(0.120 - 0.250) eV from EXO-200 [4^|42 combined with KamLAND- 

mm from the latest Planck 


ZEN 42 is represented by the horizontal dashed line, and the upper limit on 
result mi + TO 2 + TO 3 < 0.230 eV at 95% confidence level 43 is shown by vertical dashed line. 


The PMNS matrix is of the form 


UpMNS — UgrUI — 


( —i cos 9 

i cos 6 I sin 9 


1 


+ 


2 COS 9 _ sin 9 

V -\/2-\/5k 


i sin 9 cos 9 


2%5 ^/2 

K 2 sin 9 I cos 9 
2^/5 ^ x/2 / 


(3.38) 


The second column has the same form as for the GR mixing. The lepton mixing angles and 
CP phases are determined to be 

• 2/1 ^• 2/1 • 2, n 4 — 2k 

sm 6 'i 3 = ^ sm 9 , sm 9 i2 = --—^-— , 

V5 b — 2 k + cos 29 


13 











































































Figure 3: The (/3/3)oy—decay effective mass |mee| with respect the lightest neutrino mass rrimin in case IV. 
The orange and green bands denote the Scr regions for normal ordering and inverted ordering mass spectrum 
respectively. The red and purple areas are the predictions for the lepton mixing matrix in Eq. (3.381. The 
present most strict bound \mee\ < (0.120 — 0.250) eV from EXO-200 40 41 combined with KamLAND- 
ZEN 42 is represented by the horizontal dashed line, and the upper limit on rrimin from the latest Planck 
result TOi +1712 +1713 < 0.230 eV at 95% confidence level 43 is shown by vertical dashed line. Note that the 
correlation between sin^ di 2 and sin0i3 is the same as that of case III and can be found in Fig. 


sin^ 6*23 = - , |sin5(yp| = l, sin 021 = sin 031 = 0 . (3.39) 

We see that both 6*23 and 6 cp are maximal and the two Majorana CP phases 0^21 and a^i 
are trivial. Similar to case III, the relation sin^ 6*12 cos^ 6*13 = (3 — k)/ 5 is satisfied as well. 
The best fitting results for the three mixing angles are listed in Table The predictions for 
the (/9/?)o!..—decay effective mass \mee\ are shown in Fig. 

q q ^ _ yT^ST'^S n — yST'^ST^S 

0*0 Lt / — ^3 ^ 1 / — ^2 

In the charged lepton sector, the remnant CP transformation H\jp is determined by the 
consistency condition 

Xi,p:{T^ST^S)Xr,^ = Prig'), g' e . (3.40) 

We find that there are 6 possible solutions for X;r, i.e., 

H^P = {PriST^),PriST^S),PriT^),PriT^S),PriT^ST^ST^),PriT^ST^ST^S)} . (3.41) 

The charged lepton mass matrix should respect both the remnant family symmetry 
and the remnant CP symmetry H^p : 

piiT^ST^S)mlmipsiT^ST^S) = mjmi, X/^mlmiXis = (mjmi)*, X 13 E H^cp • (3-42) 

Notice that the three residual CP transformations Xir =pr{ST^), pr{T^S), pr{T^ST"^ST^S) 
lead to degenerate charged lepton masses since both p^iST^) and priT^S) are not symmetric. 
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For the remaining ones Xir =pr{ST^S), pr{T^), pr{T^ST"^ST^), the hermitian combination 
m\mi is constrained to take the following form 


m\mi 


/ a 2 + a/2(2k — 3)c) e 2Khe ^ 

2 + \/2{2k — 3)c) a + + (8k — 14)c 2(k — l)ce^ 


V 


2 Kbe^ 


,(3.43) 


2 {k — l)ce 5 * Q — ^(j) ^/ 2 c) j 

where a, b and c are real parameters. It can be diagonalized by the unitary matrix 


U, = 


( 


\ 


7-4k^-^ 
15 ^ ® 


2^ 3^ 

15 ^ ® 


\ 

15 ® ® 


Ui-dhrh* 111 + 


2 ^/l^ 

15 


i(l + 


15 

7-4k 

15 




1 ( 1 - 


7-4k 

15 


/ 


(3.44) 


with Ujm\miUi = diag(mg, m^), where the charged lepton masses are 

rUg = a — 4(k — l)c, = a — \/ 6(2 -|- k) b — ^8 — 5k + ■\/3(47 — 29k) j c , 
ml = a+ ^6(2 + k) 6 + (^5k - 8 + ^3(47 - 29k)) c . (3.45) 

The symmetry group x Hcp is broken into x H'^p in the neutrino sector. By 

solving the restricted consistency equation of Eq. (2.12a), we find 


H'^p={p,{T^),p,{TST),p,{T^ST^ST^S),p,{{ST^fS)] . 


(3.46) 


The neutrino mass matrix preserving the remnant family symmetry is of the 

form 


10 0\ /2 0 

rriy = a \ 0 0 1 ) +(3 
0 1 0 


47ri 

0 3e-— -1 


+7 


/ Sttz _ 37ri 

0 e 5 e 5 

Stti /— ■Ki 

e 5 y/ 2 e 5 0 


yO —1 3e 5 


/ ,— 27ri 27ri 

2^/2 e~~ e~ 


+6 




j - TTl 

0 \/2e 5 


27ri J — TTi ,— 

e 5 v^e 5 —\/2 
e 5 — 1/2 \/2e 5 


(3.47) 

where parameters a, /3, 7 and 5 are generally complex, and they are further constrained to 
be either real or imaginary by CP symmetry. It is convenient to hrstly perform a constant 
unitary transformation Gcpp and yield 


m!„ = UnT,pmM, 


GRpHi'u'-JGRP 

a + 2f3 — ■\/2(l -|- k) 7 


0 


0 

-a + 4/3 — \/27 0 

0 a + 2/3+ V 2 ( 2 -k) 7 , 


(3.48) 


where 


Ugrp — 




0 


- \ 

y V5 


-1 _ 3772 

72^ 

1 


K 

V5^ ^ 72 ^ ^ 


2772 


K — 


\/2V5k 
1 


e 5 


\/2V5f 


_ 2772 

e 5 


/ 


(3.49) 


Next we discuss the constraints of the residual CP symmetry on the neutrino mass matrix 


m^. 
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(V) X,, = pr{T^),pr{T^ST^ST^S) 

In this case, a, (3, 7 and 6 are determined to be real. Then neutrino mass matrix m'^ is a real 
symmetric matrix, and it can be diagonalized by a rotation matrix in the (2,3) sector. 


cos 0 0 — sin 6 
U[,= \ Q I 0 
sin 9 0 cos 9 


with 


tan 29 = 25/7 . 

The three light neutrino masses are given by 


(3.50) 


(3.51) 


1 

m, = - 


2a + 4/3 — a /27 ~ 


7107 


m 2 = 


m 3 = 


—a + 4/3 — a /27 


2a + 4/3 — a /27 + 


cos 20 


7107 


cos 29 


(3.52) 


The lepton mixing matrix is of the form 

/ cos 0 + sin 0 

UpMNS = U}UgrpUI = —^ \ COS0 - e"? sin0 

Ve¥cos0 + e¥sin0 


1 cos 0 — sin 0 

1 COS0 —sin0 

1 COS0 + sin0 


(3.53) 


We see that the second column of the PMNS matrix is (1,1,1)^/\/3, which frequently 
appears in discrete flavor symmetry models. The leptonic mixing parameters read a^ 


sin^^is = -(1 


sin 20 ), sin^ 012 = 


2 + sin 20 ’ 


sin^ 023 = 7 , 


Isin^cpl = 1, sina2i = sinasi = 0. 


(3.64) 


Both Dirac CP phase and 023 are maximal while Majorana CP phases are conserved in this 
case. In common with all trimaximal mixings, 0 i 2 and 0 i 3 are related with each other by 

3 sin^ 012 cos^ 013 = 1. (3.55) 

The measured 3cr range 0.0176 < sin^ 0 i 3 < 0.0295 gives rise to 0.339 < sin^ 0 i 2 < 0.343 
which can be directly tested by JUNO in near future 44 . The correlation between 0 i 2 and 


013 and the predictions for the (/3/3)ojy—decay are displayed in Fig. All the three mixing 
angles can agree within 3cr with the experimental data for certain values of 0. The best 
htting results are listed in Table [T| and the minimum values of the functions are 3.987 
and 7.480 for 10 and NO, respectively. 

(VI) X,, = p,{TST),p,{{ST^fS) 

The requirement of real a, /3, 7 and pure imaginary 5 follows immediately from the remnant 
CP invariant condition. In the same way as previous cases, the PMNS mixing matrix is 
found to be 


UpMNS — 


bni „ 27ri , „ 27ri „ 

e 6 cos 0 + 63 sm 0 1 es cos 0 
6 ^ cos 9 — 6 ^ sin0 
sin 9 — i cos 0 


sin0 

. 4:7ri „ Ttti , „ 

lea cos 0 + 66 sm 0 


(3.56) 


cos 9 + i sin 0 


^For cos 20 = 0, we have sin 0 i 3 = 0 or cos 0 i 2 = 0 so that 6cp cannot be determined uniquely. 
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Figure 4: Results for sin^ 6*12 and sin 6*13 (left panel) and the allowed values of the effective mass |mee| (right 
panel) in case V. On the left panel, the best fitting value is labelled with a red pentagram, and the points for 
9 = 0, 7 r /6 and 27r/3 are marked with a cross to guide the eye. The Icr and Scr ranges of the mixing angles 
are taken from Ref. |^. On the right panel, the orange and green bands denote the 3 (t regions for normal 
ordering and inverted ordering mass spectrum respectively. The red and purple areas are the predictions 
for the lepton mixing matrix in Eq. (3.53). The present most strict bound \mee\ < (0.120 — 0.250) eV from 


EXO-200 40,41 combined with KamLAND-ZEN 42 is represented by the horizontal dashed line, and the 

is 


upper limit on rrimin from the latest Planck result mi + m 2 
shown by vertical dashed line. 


- m 3 < 0.230 eV at 95% confidence level 43 


The expressions for the lepton mixing parameters are as follows, 


, n„ 1 \/3sin26' 

sm^ 013 =- 

3 6 


sin^ 023 = 


2 + a/3 sin 20 
4 + -v/S sin 20 ’ 


sin 021 


2 sin 20 + -\/3 
2 + a/3 sin 20 


Sil/ 012 = 
Isin^cpl = 
Isinogi 


2 

4 + \/3 sin 20 ’ 

8 cos 20 + a/3 sin 40 
2(2 + y3sin20)\/4-2\/3sin 20 

4a/3cos20 

5 + 3 cos 40 


(3.57) 


where = 031 — 25cp- If is remarkable that all the three CP violating phases nontrivially 
depend on the parameter 0. However, we see that in case of 0 = 7 r /4 the minimum value of 
013 is obtained with sin^ 0 i 3 |e= 7 r /4 = (2 — a/ 3)/6 ~ 0.0447 which is outside the 3 ct range j^. 
Furthermore, we note that the atmospheric angle 023 is the complementary angle of 0 i 2 or 
is equal to 0 i 2 if the second and the third rows of the PMNS matrix is interchanged. As a 
result, this mixing pattern is not compatible with experimental data and consequently we 
don’t included it in Table [H 


3.4 


tAST^ST^S,TST^) 

Cr; = A 4 , 



In the last case, the residual symmetries are assumed to be s,tst ) ^ H\jp in 

the charged lepton sector and Z/ x Tf/p in the neutrino sector. For the remnant family 
symmetry s,tst ) Pold, the mass matrix mjmi has to fulhll 


pl{ST^ST^S)m]mip3{SASTA) = m\mu pl{TST^)m\mip^{TST^) = mlmi. (3.58) 
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Then m]mi is constrained to take the form 


/ 


m\mi = 




2Khe-f 


\ 


2Kbe-"-^ 


a + 2 y/ 2 Kb + 2 {k — l)c 2 {k — l)ce 5 ' 

2 {k — l)ce^ a + 2a/2k6 + 2 {k — l)c y 


where a, b and c are real. It is diagonalized by the nnitary matrix 


/ 

V Vs 


Ui = 


\/5 K 


\/^ 

1 


_ ZttI 

-e 5 

37ri 

=e 5 


\ 

with Uim\miUi = diag(mg,m^) where 


1 _^ 

1 7 ri 


275 


27rt 

e 5 


275 


27ri 

e 5 


2a/2 6, = a + 2a/2 + 4(k — 1)c, ml = a + 2^2 n^b. 


(3.59) 


(3.60) 


(3.61) 


The mass matrix m\mi is also snbject to the constraint of the residnal CP symmetry H^p. 
It is straightforward to determine that H^p can take the valne 

Hhp= {Pr{ST^ST),p,{{ST^fS),pr{ST^),p,{T^),Pr{{T^SYT^),Pr{T^ST^),Pr{T^S), 

PriT%STy),Pr{T^ST^ST^S),pr{T^ST^),Pr{TST^S),pr{TST)}. (3.62) 

The twelve CP transformations can be classified into two categories. For X/r = pr((5'T^)^S'), 
PriT"^), Pr{T^ST"^ST^S), pr{TST), the remnant CP invariant condition Xj^mlmiXis = 
(m|mz)* is antomatically satished, and therefore no additional constraint is prodnced. Nev¬ 
ertheless, the remaining eight CP transformations Xir = pr{ST'^ST), pr{ST^), pr{{T‘^S)'^T^), 
Pr{T‘^ST‘^), priT^S), pr(T^{ST‘^y), pr{T‘^ST‘^) and priTST'^S) are not viable, as they re- 
qnire 6 = c = 0 so that the charged lepton mass spectrnm is completely degenerate with 
~ ~ ~ nentrino sector, the remnant symmetry Zf x H^p and its phe¬ 

nomenological implications have been stndied in section in The nentrino mass matrix mi, 
is fonnd to be given by Eq. (3.8), where the parameters a, f3 and 7 are real while 6 is real 


or pure imaginary depending on the residual CP transformation Xy,.. 

(VII) X^r = Pr{T^ST^ST^),p^{T^ST^ST^S) 

In this case, the neutrino mass matrix is diagonalized by the unitary matrix in Eq. (3.21). 


Combining the unitary transformation Ui in Eq. (3.60) from the charged lepton sector, we 
obtain the lepton flavor mixing matrix: 


UpMNS— 2 


K cos6 '-|-(k — 1) sin6' (k — 1) cos 0 — sin 6^ 
— 1 (k — 1) cos6*-|-fi:sin6' kcos^ — (k — 1) sin6' 
K — 1 sin 9 — K cos 6 cos 6 + n sin 9 


(3.63) 


where the parameter 9 is specified by Eq. (3.19). The lepton mixing parameters are predicted 
to be 


sin^ = 


cos 9 — Ksin 9Y 

13 = --) 


sin^ 003 = 


4/s:2 

{k^ cos 9 — sin 0)^ 

Ak'^ — (cos 0 — K sin 0 ) 2 ’ 


sin^ 012 = 


(KCOS 0 -|- sin 0 )^ 


4^2 — (cos 0 — Ksin 0)2 ’ 
sinbcp = sin 021 = sin 0:31 = 0 . 


(3.64) 


18 





















We find all the three CP violating phases 5cp, 021 and ctai are conserved, this is becanse that 
a common CP transformation pr{T^ST"^ST^S) is shared by the nentrino and charged lepton 
sectors. In addition, 623 deviates from maximal valne. After some tedions calculations, we 
hnd the following relations between the mixing angles 


4 cos^ 612 cos^ 013 = 1 + K , 

5 sin^ 023 = 3 — K + (1 + 2k) tan^ 0i3 ± 2k tan 0i3 + k — (2 + 3k) tan^ 0i3 , (3.65) 

which is plotted in Fig. For the 3a interval 0.0176 < sin^ 0i3 < 0.0295 [^, we have 0.326 < 
sin^ 012 < 0.334 and 0.454 < sin^ 623 < 0.511, which are in the experimentally favored 
ranges [^. The global minimum of the function is rather small 3.503 (1.626) for NO 
(10) neutrino mass spectrum, therefore this mixing pattern can describe the experimental 
data very well. Moreover, we note that the best fitting value of 623 is in the first octant 
with sin^ 023(^6/) = 0.480 (0.486) for NO (10) spectrum. Agreement with experimental data 
can also be achieved if the second and third rows of the PMNS matrix in Eq. (3.63) are 
exchanged. Then the atmospheric mixing angle 623 changes to 

'^(cos0 + Ksin0)^ 


sm ^23 


023 = 


K 


4k‘^ — (cos0 — fi:sin0)^ 


(3.66) 


and the predictions for the other mixing parameters remain as Eq. (3.64). The allowed region 
of sin^ 023 becomes 0.489 < sin^ 023 < 0.546 with the best fitting value sin^ 023(0b/) = 0.510 
(0.513) for NO (10) spectrum. Obviously 023(0b/) is in the second octant. Comparing with 
other mixing patterns shown in Table [T| we see that this case gives rise to the smallest 
Xmin both NO and 10. The above predictions for solar and atmospheric mixing angles 
could be tested directly in near future, since the next generation neutrino oscillation exper¬ 
iments are expected to reduce the experimental error on 0i2 and 623 to few degrees. The 
theoretical results for the (/9/3)oi/—decay effective mass \mee\ are displayed in Fig. Note 
that interchanging the second and third rows does’t matter since \mee\ is independent of 623- 
Again, the predictions for 10 neutrino spectrum are within the sensitivity of forthcoming 
experiments. 


(VIII) X,r = Pril),Pr{S) 

The neutrino mass matrix is diagonalized by the unitary transformation in Eq. (3.15). The 
PMNS matrix is found to take the following form 

^ / sin0 —mcos0 cos0-|-zfi:sin0 k —l\ 

UpMNS = - \ icosO + {k — l)sm6 (k — 1) cos0 — isin0 k , (3.67) 

y i(fi: — 1) COS0-1-fi:sin0 k cos 0— *(/? —1) sin 0 —1 J 

up to permutations of rows and columns. The lepton mixing angles and CP phases can be 
read off as 


sin^ 013 = 


Isinhcpl = 


75 


~ 0.0955, sin^ 0i2 = ^ ^ cos 20, sin^ 023 = 


VTOsin20 


79 — cos 40 


sma2i = 


2 10 
8 sin 20 


9 — cos 40 


isinog]^! = 


5 + 75 
10 

2 sin 20 


~ 0.724, 


75 + cos 20 


(3.68) 


We see that the solar mixing angle 0i2 has a lower bound given by sin^ 0i2 > (5 — \/5)/10 ~ 
0.276, and the experimental data on 0i2 can be accommodated for particular values of 0. 
Both 013 and 623 are independent of 0, and they are outside the 3 ct ranges [^. Furthermore, 
6x6 =36 possible permutations of rows and columns of this mixing pattern are considered. 
However, none of them can give rise to three mixing angles in the experimentally preferred 
3a range [^. 
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Figure 5: The results for sin^ 0 i 2 , sin^ 6*23 and sin 6*13 (the former three panels) and the allowed values of 
the effective mass \mee\ (the last panel) in case VII. The global minimum of the function is labelled with 
a red pentagram, and the points for 0 — 0, 7 r/ 6 , 7 r/ 3 , tt/2, 27r/3 and 5tt/6 are marked with a cross to guide 
the eye. The black solid lines and blue dashed lines in the upper-right and lower-left panels represent the 
two solutions for 023 shown in Eq. (3.64) and Eq. (3.66) respectively. The corresponding PMNS matrices 
are related through a exchange of the second and third rows. The Icr and 3cr ranges of the mixing angles 
are taken from Ref. [^. In the last panel, the orange and green bands denote the 3a regions for normal 
ordering and inverted ordering mass spectrum respectively. The red and purple areas are the predictions 
for the lepton mixing matrix in Eq. (3.63). The present most strict bound \mee\ < (0.120 — 0.250) eV from 


EXO-200 [4I)||44] combined with KamLAND-ZEN 42 is are represented by the horizontal dashed line, while 
the upper limit on rn^m from the latest Planck result mi +m 2 +m 3 < 0.230 eV at 95% confidence level 43 
is shown by dashed line. 


4 Model building 


In previous section, we have performed a model-independent analysis of the lepton mixing 
patterns which can be derived from xi Hcp- As summarized in Table [T| we hnd hve 
new mixing patterns which are compatible with current experimental data. In this section, 
we shall construct a concrete model with both A 5 family symmetry and generalized CP 


symmetry, the symmetry breaking patterns studied in section 3.2 are implemented, and 


therefore the lepton flavor mixings given by Eqs. (3.32, 3.38) in case III and case IV are 


realized. Note that it would be also interesting to implement other cases such as case VII in 
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Table 2: The matter fields, flavon fields, driving fields and their transformation properties under the family 
symmetry A 5 x Z 3 x Z 4 x and C/(l)i^, where the phase wa = and we = e^. 


a model. In the present model, both the three generations of left-handed lepton doublets I 
and the three generations of right-handed neutrinos are assigned to transform as triplet 
3, while the right-haned charged leptons and are all invariant under A^. In discrete 

flavor symmetry model building, either cyclic Zn or continuous U{1) symmetry is frequently 
introduced to eliminate unwanted operators, to ensure the required vacuum alignment and 
to reproduce the observed charged lepton mass hierarchies. The auxiliary symmetry is taken 
to be Z3 X Z4 X Zq in this model. The A^ family symmetry and CP symmetry are broken by 
some flavons in a proper manner. All the flavon helds are standard model gauge singlets. As 
anticipated, we formulate our model in the framework of supersymmetry (SUSY). A U{1)r 
symmetry related to i?—parity and the presence of driving helds in the havon superpotential 
are common features of supersymmetric formulations. The held content of the model and 
their classihcation under the symmetry are listed in Table In the following, we hrst discuss 
the vacuum alignment of the model, then specify the structure of the model at leading order 
and next-to-leading order. As we shall show, the lepton mixing is exactly the GR at LO, and 
a non-vanishing value of the reactor mixing angle ^13 is generated by higher order corrections. 
Consequently 6^13 is naturally of the correct order in our model. 

4.1 Vacuum alignment 

We utilize the standard supersymmetric driving held mechanism to solve the vacuum 
alignment problem. A global U{l)ii continuous symmetry is assumed in this approach, and 
the usual R—parity is a discrete group of this U{1)r. The matter helds have R—charge 
equal to one, both havon helds and Higgs are chargeless and the driving helds carry two 
units of R—charge. At LO the most general driving superpotential Wd invariant under 
A 5 X Z 3 X Z 4 X Zq with R = 2 can be written as 

= (4.1) 


with 


+ /2(0°(</30)4)i + /3(0°(V^^)4)1 + M^(VzV)l + /4(^°(w)5)l , (4.2) 
Wd = + g2^°{xx)i + 93^°{pp)i + g4^{x°x)i + ^5(x°(xA)3)i 

+5'6(p°(p2k)3)i -|- Ma(A°A)i -|- gj{A^{xx)5)i + 5'8(A°(pp)5)i, (4.3) 


where (.. .)r denotes a contraction into the A^ irreducible representation R according to 
the Clebsch-Gordan coefficients listed in Appendix]^ Notice that all the couplings fi{i = 
1,..., 4), gfj (z = 1,..., 8) and the mass parameters M^, M^, Ma are real, since the theory 
is invariant under the generalized CP dehned in Eq. (2.9). In the SUSY limit, the vacuum 
alignment is achieved via the requirement of vanishing F—terms of the driving helds. In the 
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charged lepton sector, the equations for the vanishing of the derivatives of with respect 
to each component of the driving helds are: 


da^ 

dw\ 

d(l)\ 

dwh 

W 2 

dwh 

M 

dwh 

M 

dwh 


dilji 


= + 2(P2(^3) = 0, 

= /2(</2203 + \/2(p30i) - /3(2\/2(pi'05 + </?2V’4 “ \/6(p3'0i) = 0, 

= ~/2(\/2</3i 03 + </3202) + /3(V^</3l'04 + ~ = 0, 

= -/2(\/2<Pi02 + ^3h) - /3(\/2(pi^/’3 - 2ip2'lp2 + ^<^3^^^) = 0, 

= /2(\/29?201 + ‘^3<t>2) + h{2,\/2ipi'lp2 - \/6(/92'01 + </93l/>3) = 0, 


= + 2 fi{(pl - IP2(P3) = 0, 


dwl 


d^jl 

dw'd 

d'lps 

9w^d 

dw'd 

d'lpl, 


- 2v/3/4(pi(p3 = 0, 

= M^ipi + '/Qf4^pl = 0, 


= M^V’3 + = 0, 

= M^'iIj2 - 2a/3/4(Pi(P2 = 0 , 


We hnd one solution to those equations, 


(‘d) = 






( 0 ) = 






(^) = 


/0\ 

0 

1 

0 

W 


: 


( 4 . 4 ) 


( 4 . 5 ) 


up to family symmetry transformations, where the vacuum expectation values (VEVs) 
v^p, V(j, and are related by 




375 / 3/4 2 
M^/2 


75/4 


(4.6) 


with undetermined. A common order of magnitude for the VEVs (scaled by the cutoff A) 
is expected. In order to generate the mass hierarchies among the charged lepton, we assume 


X 


rsj 


- r\j - 

A A 


0 {X 


(4.7) 


where Ac — 0.23 is the Cabibbo angle (^. In the neutrino sector, the minimization equations 
for the vacuum are 


dxv^ 

+ giCp + g2{xl + 2 x 2 X 3 ) + 93{pI + ‘2p2p3) = 0, 

= 5'7Xi - 5'5(2xiAi - 73 X 2 A 5 - V3x3^2) = 0, 
dXi 
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dw'i 

dw’i 

dw'i 

M 

dw'i 

M 

M 

aA? 


8A« 

SA« 

Mf 


= 9iiXz + 5'5(\/3xiA5 - VQx2^4 + X 3 A 1 ) = 0 , 

= 9i^X2 + 5'5(\/3xiA2 + X 2 A 1 - ^ 6 x 3 ^ 3 ) = 0, 

= ge{V3piAi + P 2 A 4 + P 3 A 3 ) = 0 , 

= 5'6(piA5 — V 2 P 2 A 3 — V 2 P 3 A 2 ) = 0, 

= 9eiPi^2 — V 2 P 2 A 5 — V^p^A^) = 0, 

= MaAi + 2 g ' i{xi ~ X2X3) + 2(y'8(Pi “ P2P3) = 0, 
= MaAs - 2 ^ 85 ( 7 X 1 X 3 + VQPspI = 0 , 

= MaA4 + 1 / 657 X 3 - 2a/358PiP3 = 0, 

= MaA3 + 1/657X2 - 21/358P1P2 = 0, 

= MaA 2 - 2 ^ 857 X 1 X 2 + V 658 P 3 = 0 . 


A solution to those equations with each flavon acquiring non-zero VEV is given by 


/ 


A1' 

(0 = (C) = v^, (x) = I 1 I (P) = 



KVi 
Vl 

-(1 + K.)vi 
-(1 + k)vi 
V l 


\ 


These VEVs are related through 

10 (^- 8)5557 2 
= -TV- v: 


94 MA 


-A 


2(« - 8) [(5258 + 9397) 9 iMA - 5555758M5] 


9i9a9sMa 


2 (2 1^)97 2 

V = - V. 

98 


X’ 


Ui = 


1 / 30(2 - k )57 

Ma ' 


X ’ 


(4,8) 


(4.9) 


(4.10) 


where is undetermined. It is easy to check that the VEVs of ( and A break the A 5 family 
symmetry down to ^ while the subgroup is preserved by vacuum of 


X and p. Furthermore, Eq. (4.10) implies that and Ui have the same phase 

up to TT, since all couplings are real. In our model, the GR mixing is reproduced exactly 
and a non-zero reactor mixing angle ^13 is generated after subleading order contributions are 
included. In order to obtain the correct size of ^ 13 , we could choose 


U 5 uc Vp Vl 

- rsj - r\j - rsj - rsj - r\j I y I A ^ 

A A A A A ^ 


(4.11) 


4.2 Leading order results 

The charged lepton mass terms, which are invariant under the imposed family symmetry 
A5 X Z 3 X Z4 X Zg, can be written as 

wi = ^T‘'{lp)ihd + ^/i/^(0/’)3)ihd ^/i//(^/’^)3)ihd ^e''(/(p)i(00)ihd 
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( 4 . 12 ) 


+ ^e^{l^)i{.'ip'ip)ihd + ^e^{(l^)s{'ip'ip)3,)ihd + ^e^((V)5(V''0)5i)i/id 

+ ^e'^((/99)5(V’V’)52)l^rf + ••• ) 

where dots stand for higher dimensional operators which will be discnssed later. Note that 
all conplings here are real dne to the generalized CP symmetry. After the electroweak and 
flavor symmetries breaking by the VEVs shown in Eq. (4.5), we obtain a diagonal charged 
lepton mass matrix, and the three charged lepton masses are 

rrie = \/2 


q ^ / 

3|/e2 A Q A (l/es 


= a/2 


2 /mi 


A3 




A3 


A3 


Vd, 


A2 




v^ 

Vd, 

m^ = 



Vd, 


(4.13) 


We see that the realistic mass hierarchies me : : m,- ~ : 1 is generated for the 

order of magnitnde of the flavon VEVs in Eq. (4.7). Fnrthermore, as both mi and psiT) are 
diagonal, obvionsly we have pl(T)m]mip 3 (T) = m\mi^ i.e., the residnal flavor symmetry of 
m\mi is Z/. Next let’s discnss the nentrino sector. Nentrino masses are generated by type 
I see-saw mechanism in this work. The LO snperpotential for nentrino masses is 




( 4 . 14 ) 


where the conpling constants |/i, 1/2 and the mass M are enforced to be real by the generalized 
CP symmetry. The Dirac mass matrix is obtained from the first two terms in Eq. (4.14) and 
it is given by 


( 4 . 15 ) 


lOoA / 

^-272k 

-3 

-3 

001 Wn + 6 j 

-3 

-372(k + 1) 

a/2k 

0 10/ \ 

. -3 

y/2K 

-37^(fi; + 1) 


where Vu = (hu), and the parameters a, b are 

a = yiir, b = y2 


(4.16) 


A’ “ ""ysA' 

The common phase of a and b can be absorbed by field redefinition, conseqnently both a and 
b can considered as real. The last term of Eq. (4.14) leads to the Majorana mass matrix: 

( 4 . 17 ) 


1 0 0 

mM = M \ 0 0 1 


0 1 0 

Therefore the three right-handed nentrinos are completely degenerate with mass eqnal to 
M. The light nentrino mass matrix is then given by the see-saw relation: 


m^, = —m'^mjjmf) = a 


/lOO^ 

l+7f 

0 0 1 

[oloj 

1 V2 \ 


- 2 V 2 3 3 \ /2 0 0 

3 072 +7 03 - 1 ], (4.18) 

3 720 / \0-13 


where 


a 


— - [7 + 40(1 + k)7] 
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(4.19) 


/3 = 2 \V2ab - {3 + 4:K)b^ 


K 

M’ 


7 = 


272(1 + ^)06+ (1 + 8 k )6^ 


K 

M 


We find that the neutrino mass matrix rrty in Eq. (4.18) is of the same form as the general 
mass matrix in Eq. (3.27) with 6 = 0. Therefore rrii, is exactly diagonalized by the GR 
mixing pattern, i.e., 

Uqr^uUgr = diag(mi, m2, m3), (4.20) 

where the phase matrix Ki^ which encodes the CP parity of the neutrino state, has been 
omitted. The mass eigenvalues mi 2,3 are 

,2 


mi = 

m2 = 
m3 = 


- 272(3 - R)ab + 10(2 - k)}/ 


M’ 


— 10y/2Rab + 50(1 + k) 6^ 


M’ 


7 + 272(3 + 4706 + 10(5 + 876^ 


M 


(4.21) 


Since the charged lepton mass matrix is diagonal in LO, the lepton mixing is exactly the GR 
mixing pattern. Here the reason why the GR mixing is produced is because that the flavor 
symmetry is broken to K\ ’ subgroup by the VEVs of ^ and A. Furthermore, 

we see that three neutrino masses mi^2,3 only depend on two real parameters a and 6 which 
can be hxed by the measured values of the mass-squared difference 6m^ = m\ — m\ and 
Am^ = rn\ — (m^ + m|). For the best htting values 5m? = 7.54 x 10“®eV^ and Am^ = 2.43 x 
10“^eV^ we hnd the neutrino mass spectrum can only be NO, and the absolute values 
of the neutrino masses are mi = 4.81 x 10“^eV, m2 = 8.70 x 10“^eV and m3 = 0.0497eV. 


4.3 Next-to-leading-order corrections 


At LO our model gives rise to the GR mixing pattern Ugr which predicts a vanishing 
reactor mixing angle (^13 = 0°). Hence substantial next-to-leading-order corrections are 
needed to bring the model to agree with the experimental data on 6^13. We will demonstrate 
in the following that a non-zero 6^13 can be obtained after the NLO contributions are included. 
Moreover, the LO remnant symmetry K\ ’ ^ of neutrino sector is further broken d own 

to mixing patterns of case HI and case IV discussed in section 

are realized. Firstly we consider the corrections to the flavon superpotential in Eq. (^ 


3.2 


. 2 ) 


which determines the vacuum alignment of the charged lepton sector. The symmetry allowed 
NLO operators are of the following form 

+ = ((/v)5(w)5)i/A + , (4.22) 


where all possible contractions should be considered, and all dimensionless coupling 
constants are omitted with = {cr°,'0°}, T; = {0,'0}, = {7^} and = {C,x}- 

Note that 5w\ is suppressed by with respect to the LO superpotential in Eq. (0. 
The NLO vacuum conhguration is determined by searching for the zeros of the F—terms of 
+ 5w^^ with respect to the driving helds ct°, (jP and AVe hnd that the NLO vacuum of 
+, 0 and 0 are given by 


(+) 


/ egA^ \ 


oAn 

/ \ 


€ 7 X 1 

1 

(0) = 1 + e^Xl v^, 

(0) = 

1 + cgA^ 

oA^y 

\ J 


CgA^ 

\ eioAc ) 


(4.23) 
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where e* (i = 1,..., 10) are general complex numbers with absolute values of order one. The 
higher dimensional operators contributing to the charged lepton masses are: 




(4.24) 


The charged lepton mass matrix can be obtained by inserting the NLO VEVs of Eq. (4.23) 
into the LO mass terms plus the contribution of 6wi evaluated with the LO VEVs of Eq. (4.5). 
We find that the NLO charged lepton mass matrix is of the following form: 


mi ~ 


/ me A^me A^meX 
\ X^mr X^mr mr ) 


(4.25) 


Therefore the contributions of charged lepton sector to the lepton mixing angles is of order 
Ae and can be neglected. 

We proceed to discuss the subleading corrections in the neutrino sector. The higher order 
corrections to the flavon superpotential of C, y, p and A read 

= f c^(x\)i + ^-fax\xxh)i + ^(x\)i(xx)i + ^-f{{x\Uxxh)i 

+ ^((x°x)5(xx)5)i + ^(x°x)i(pp)i + ^((x°x)5(pp)5)i 

+ ^((p°p)5(xx)5)i + p)h{pp)^)i 1 (4.26) 


where all couplings Pi (z = 9,..., 17) are real because of the generalized CP symmetry. The 
resulting contributions to the F—terms of the driving helds a^, and are suppressed 

by (T)/A Ac (4/ = {(C, X, p}) compared to the contribution from the LO terms in Eq. ( |4.3[ ). 
Hence they induce shifts in the VEVs of (, y, p and A at relative order Ac with respect 
to the LO results. After some straightforward algebra, the new vacuum configuration can 
be written as 


where 


(0 =v^ + 6v^, (C) =v^ + Sv^, (x) = 

/ 


-y/2K{Vp + 6Vp) 

(p) =1 VpA 5vp 

Vp + 6vp 


(A) = 


K A. 


3 {-KVi + (1 + 2k) 6va) 
vi + 6va 

(1 + k) ui + 2k6va 
(1 + k) ui + 2kSva 
vi + 6va 


\ 


5vf = 


6vc = 


5vp = 


Xi + ^(5X2 
5'45'6A 

2{k- 2)giQMAvl + 2gnMAvl 
Ag^pfiAvp 


gsM^Xi + {g^gsM^ — g^g4:MA)X2 


2fi'i5'4fi'65'8Anf 


5va = 


2a/6((k - 1)5-16^^ - ngnvl 

PeA 


(4.27) 


(4.28) 


with 

Ai = ge + 2(3 — K){gii + ^gi^)v'^ + 2\/hK{gi4^ + gi^)v^^ , 

X 2 = 2{k - 3)^716^2 + 2V5t^gi,vl . (4.29) 


26 













Obviously the vacuum of x is kept intact, (p) acquires 0{Xc) corrections in the same direction, 
while the alignment of A is tilted. Moreover, from the relations in Eq. (4.10), we see that the 
shifts 6v^, 6v(^, 6vp and 6va carry the same phase as v^, Vp and vi up to tt, respectively. 

The light neutrino mass matrix receives corrections from both the modihed vacuum and 
the higher dimensional operators in the superpotential Wy. It is easy to check that the 
NLO corrections to the Majorana mass terms are suppressed by 1/A^ which can be safely 
neglected. The subleading operators contributing to the neutrino Dirac masses are as follows 

+ + ^{li^^)i{pp)ihu + ^{{li^'')5{pp)5)ihu ■ (4.30) 

As a consequence, the corrected Dirac mass matrix becomes 


= a 


V -3 -V2 



3 -3 

0 -V2 


Vu + C 


0 



0 0 
3^2 -^/2 

- V2 3^2 


(4.31) 


where the four parameters a, b, c and d are 


a = yi 

b = 


Vf + 6v£ n K 

- X-^ ^ A2 + 275^1/8^, 

y2 vi + 5 va 


c = 


73 A 
1/2 (1 + k)vi- 2k6va 
73 A 


2 2 

-272(k-1)//7^ + 72//9^, 




d = yA 


A2 


(4.32) 


Notice that the three parameters a, b and c have the same phase with up to tt, while the 
phase difference between d and is 0, tt or ±| depending on the product PiMa [((/ 25'8 + 
5 ' 35 ' 7 ) 5 ' 45 ' 8 ^a — 5(/4P5(/7(/|Mg] being positive or negative. Since the phase of can be fac¬ 
torized out as an overall phase of the neutrino mass matrix rriiy, the VEV can be taken 
to be real without loss of generality. As a result, a, b and c are all real and the param¬ 
eter d is real for QiMa [{9298 + 9397)9198Ma - ^ 9 A 9 b 979 lM^ < 0 or pure imaginary for 
9iMa [{9298 + 9397)9198Ma — 5(/4(/5(/7p|M^] > 0. In addition, we see that d are suppressed 
by Ac with respect to a, b and c, i.e., 

a 1) ^ c ^ 0{Xc), d ~ 0{X^). (4.33) 

Utilizing the see-saw formula, we find the light neutrino mass matrix rrii, is of the same form 
as Eq. (3.27) with 


a 


= - [3a^ + 24(27 + bc + 27) - 4(3 - /y)7] 


3M’ 


/3 = 


Qb{V2a + b + 4c) — 2 {k — 1)7 


3M’ 


7 = - 


672ac + 3(36 + 2c)(6 - 2c) + 7 


3M’ 


27 



















6^-2,d[b + 2{K-l)c]'^ 


( 4 . 34 ) 


Note that the term proportional to 5 spoils the LO GR mixing, and it is of relative order 
Ac compared with a, /3 and 7 since it is induced by the NLO corrections. Therefore the 
correct size of the reactor mixing angle ^13 can be naturally achieved in our model. After 
extracting the overall phase of the parameters a, (3 and 7 are real while 5 is real or pure 
imaginary. In the case of giM^ [{9298 + 9397)9198^^ — < 0, 5 is real such that 

the neutrino mass matrix rriy has the most general form compatible with the preservation of 
the remnant symmetry x H^p with H^p = {pj.{l), Pr{T^ST‘^ST^)}. This is the 

case III investigated in the model independent analysis of section |3.2[ The lepton mixing 
matrix Upmns and the corresponding preditions for the lepton mixing parameters are given 
by Eq. (3.32) and Eq. (3.33) respectively. There is no CP violation in this case. 

In the case of giM^ [{9298 + 9397 ) 9498 ^^ — ^94959798^^] > 0; the parameter 6 becomes 
imaginary. The origin symmetry A 5 x Hcp is broken down to x H^p with H^p = 

{Pr{S), Pr{T^ST"^ST^S)} in the neutrino sector. The neutrino mass matrix rrii, has the same 
form as that of case IV discussed in section 3.2 Both atmospheric mixing angle and Dirac 


CP phase are predicted to be maximal while Majorana CP phases are conserved, as shown 
in Eq. (3.39). In short, our model reproduces the GR mixing at LO, and realistic value of 
013 is obtained after higher order contributions are taken into account. Depending on the 
overall sign of the product giM/^ [{9298 + 9397 ) 9498 ^^ — , either case III or case 

IV can be realized. 


5 Conclusions 


Combining a discrete flavor symmetry with a CP symmetry is a very promising approach 
of predicting both lepton mixing angles and CP phases. In this work we have performed a 
comprehensive analysis of the A 5 family symmetry and CP symmetry. Since the inverse of 
each conjugacy class of A 5 is equal to itself, all the inner automorphisms of A 5 are class- 
inverting while the unique nontrivial outer automorphism of is not. As a result, the 
physical CP transformations are dehned by the inner automorphisms of A 5 . In our working 
basis, the CP transformations are found to be of the same form as the flavor symmetry 
transformations. 

Assuming neutrinos are Majorana particles, we have analyzed the possible symmetry 
breaking patterns of A 5 x Hcp and the corresponding predictions for the PMNS matrix 
as well as the lepton mixing parameters in a model independent way. We hnd hve phe¬ 
nomenologically interesting mixing patterns summarized in Table and one column of the 


PMNS matrix is hxed to be 


\T t ^ tJ- J_ J_ 

! 1 \\ \ 2 ./¥.i \ 2 ./¥.) ’ C/s’ ,/s’ ,/a 


V V ’ ^x/3’ V3’ x/3^ 

or (|, —where k = (1 -|- \/h)/2 is the golden ratio. All the three mixing angles are 
determined in terms of a single real parameter 0 , and their measured values can be accom¬ 
modated for certain values of 0. In particular, the Dirac CP violating phase 5cp is predicted 
to be trivial or maximal while the Majorana phases are trivial. In contrast, 5cp is quite 
weakly constrained and Majorana phases can not be predicted if CP symmetry is not con¬ 
sidered, as shown in Appendix Our theoretical predictions can be tested by forthcoming 
long-baseline neutrino oscillation experiments such as LBNE, LBNO and HyperKamiokande. 
The predicted mixing patterns would be ruled out, if signihcant deviations of 5cp from triv¬ 
ial and maximal values were detected. Furthermore, the phenomenological predictions for 
the (/d/d)o!y—decay are investigated. The present experimental bounds are saturated, and the 
effective mass \mee\ is found to be within the sensitivity of future decay experiments 















for inverted ordering nentrino mass spectrnm. 

Gnided by above model independent analysis, we construct a flavor model with both 
A 5 flavor symmetry and generalized CP symmetry. The lepton mixing is exactly the GR 
pattern at LO, the observed mass hierarchies among charged lepton are generated, and 
the three light neutrino masses effectively depend on two real parameters which can be 
hxed by the measured values of the mass-squared splittings. Therefore the neutrino mass 
spectrum can only be normal ordering and the absolute neutrino masses are predicted. 
The model is built in such a way that the GR mixing is modihed by NLO contributions 
and only the second column of GR mixing matrix is kept. A non-vanishing value of 6*13 
is generated at NLO and it is naturally of the correct order Ac in our model. In case of 
giM^ [(fi'25'8 + 92.97) - ^g^g^giglM^ < 0, Dirac CP phase 5cp is 0 or vr, consequently 
the mixing pattern of case III of general analysis in section PI is reproduced exactly. In case 
of giM/^ [{ 929 ^ + 5 ' 3 fi' 7 ) 5 ' 45 ' 8 AfA - ^gAg^giglM^] > 0, Dirac CP phase 5cp is maximal with 
5cp = ± 7 r/ 2 , the mixing pattern of case IV is generated. In other words, our model provides 


an explicit dynamical realization of the assumed symmetry breaking pattern in section 3.2 


It is interesting to implement any of the remaining cases II, V and VII in Table in a 
concrete model. Moreover, the group X', which is the double cover of A 5 , may deserve to 
be studied in a similar fashion. Since X' has doublet representations [^, quark masses and 
mixing should be easily reproduced. 
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Appendix 


A Group Theory of 

A 5 is the group of even permutations of five objects, and it has 5!/2 = 60 elements. 
Geometrically it is the symmetry group of a regular icosahedron. group can be generated 
by two generators S and T which satisfy the multiplication rules (4^ : 

52 = T® = (5T)3 = 1. (A.l) 

The 60 element of A^ group are divided into 5 conjugacy classes: 


IGi: 
15^2 : 


20 ^ 3 : 


12 ^ 5 : 

12 G^ 


ST'^ST^S, TST\ T^{STY, T^ST\ {T^S)^T^S, ST^ST, S, T^ST'^ST^, 
T^ST^ST^S, T^ST\ T^ST^ST^S, TST^S, ST^ST^S, T^ST, {T^SfT^ 

ST, TS, ST\ T^S, TST\ T^ST\ T‘^ST\ T^ST, T^ST^, T^ST^, TST^S, T^ST^S, 
T^ST^S, ST^ST^ ST^ST, ST^ST^, {T^SfT\ T\T‘^Sf, {ST^)‘^S, {ST‘^fT‘^ 

T, T^ ST^, T^S, ST^, T^S, STS, TST, TST^T^ST, T^ST\ T^ST^ 

T 2 , j . 3 ^ {ST^, {T^Sf, {ST^y, {T^Sf, {T^SfT^, 

T\ST^)\T^ST‘^ST\ T^ST'^ST^, (A.2) 


where nCk denotes a class with n elements which have order k. The group structure of A 5 
has been elaborately analyzed in Ref. 47 . Following the convention of Ref. 47 , we find that 


A 5 group has thirty-six abelian subgroups in total: fifteen Z 2 subgroups, ten Z 3 subgroups, 
five 7^4 subgroups and six subgroups. In terms of the generators S and T, the concrete 
forms of these abelian subgroups are as follows: 


• Z 2 subgroups 

Zf^ST^s ^ ST‘^ST^S}, Z^^^" = { 1 , TST"^}, = { 1 , T^ST^}, 

Z^^ST^ = { 1 , T^ST^}, ^ {T^SfT^S}, = {1, ST^ST}, 

Z! = { 1 , 5}, = {1, T^ST^^ST^}, = {l,T^ST‘^ST^S}, 

Z^'ST^ = { 1 , T^ST^}, = {1, T^ST'^ST^S}, = { 1 , TST^S}, 

^ST^ST^s ^ ST^ST'^S}, = {1, T^ST}, zf = {1, {T'^SfTA- 


All the above fifteen Z 2 subgroups are conjugate to each other. 

• Z 3 subgroups 

^T^ST^S ^ t^ST^S, ST^ST^}, zp^"^ = { 1 , TST^S, {ST^T^}, 
ZT^ST ^ T^ST^}, Z^^ = {1, ST, T^S}, 

= { 1 , {T^SfT^, [ST^fS], Zj^^" = {1, TST\ T^ST^}, 

= {1, T^FT^, T^ST^}, Zl^ = {1, TS, ST^}, 

ZST^ST ^ ST^ST, T‘^{T‘^Sf}, = { 1 , ST‘^ST\ T'^ST^S}. 


The ten Z^ subgroups are related with each other by group conjugation. 
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• K 4 subgroups 


KfT^ST^S,TST^) ^ ^ST^ST^S ^ ^TST^ ^ {1, ST^ST^S,TST^^T^iST^}, 
Kf^ST^ST^ST) ^ ^t25t3 ^ ^ST^ST ^ s^i^T^ST^, {T^SfT^S, ST^ST}, 
KfT^ST^ST^) ^ ^ ^T^ST^ST^ ^ {1, S,T^ST^ST^,T^ST^ST^S], 

J^{T^ST\TST^S) ^ ^ ^TST^s ^ {I^ T^st"^, T^ST"^ST^S, TST"^S}, 

KfT^ST^S,T*ST) ^ ^ST^ST^S ^ ^T*ST ^ [i^ ST^ST‘^S,T^ST, {T'^SfT^}. 

All the five K 4 subgroups are conjugate as well. 

• Z 5 subgroups 

= {1, STS, ST^S, ST^S, TST], = {1, ST^, T^S, {ST^f, {T^Sf}, 

ZT^ST ^ i^i^T^ST, T^ST^, T3(5T2)2,T^AT^AT^}, zJ = {1, T, T^, T^}, 

ZTST^ = {l^TST‘^,T^ST^,{T‘^SfT^,T^ST^ST^}, = [I, ST'^,T^S,{ST‘^f ,{T^Sf}. 


All the six Z 5 subgroups are related to each other under group conjugation. 


Here the superscript of a subgroup denotes its generator (generators). The A 5 group has five 
irreducible representations: one singlet representation 1 , two three-dimensional representa¬ 
tions 3 and 3', one four-dimensional representation 4 and one five-dimensional representation 
5. In the present work, we choose the same basis as that of Ref. 47 . The explicit forms of 


the generators S and T in the five irreducible representations are as follows 


1 : 

3 : 

3' : 

4 : 


5 : 5 


5 = 1, 

1 - %/2 - y 2 ^ 

5 = ^1 -V 2 -K K-1 

-a/2 k — 1 — K 



/ 

-1 y2 

V2 \ 



5 = 

A 

y2 i-K 





a/2 k 

1 — kJ 




/ 1 

K — 1 

K — 

1\ 


Q — -L 

K — 

K 

1 -1 

1 

1 K 

— 1 K — 

1 



V -1 

K K 

- 1 1 



/-I ye 

ye 

ye 


ye \ 

a/6 (k 

-1)^ 

— 2k 

2(k- 1) 



ye - 

- 2k 


{K-lf 

2( 

K — 1) 

ye 2 {k-i) 

1 

to 



— 2k 

vye 

K? 

2{k-1) 

— 2k 

1 

to 


T = l, 

/I 0 0\ 

T — j 0 cus 0 I , 

Vo 0 ulj 


I 0 0\ 

T = 0 a;2 0 , 

Vo 0 ulj 


T 


T 


( CJ 5 0 0 0 ^ 

0 cn2 0 0 

0 0 cnf 0 

Vo 0 0 ujI) 

/I 0 0 0 0 \ 

0 cus 0 0 0 

0 0 0 0 
0 0 0 cnf 0 

yo 0 0 0 LolJ 


(A.3) 


2772 nn 

where U 5 = . The character table of A 5 group is reported in Table We can straight¬ 

forwardly obtain the Kronecker products between various representations: 


1(8)R = R(8)1 = R, 3®3 = 1©3©5, 3'©3' = 1©3'©5, 3x3'= 4© 5, 

3®4 = 3'©4©5, 3'©4 = 3©4©5, 3©5 = 3©3'©4©5, 

3'©5 = 3© 3'©4© 5, 4©4 = 1©3©3'©4©5, 4©5 = 3©3'©4©5i©52, 

5 © 5 = 1 © 3 © 3' © 4i © 42 © 5i © 52. (A.4) 
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Conj 

ugacy Classes 

IGi 

ISC's 

2 OC 3 

I 2 C 5 

12 C' 

1 

1 

1 

1 

1 

1 

3 

3 

-1 

0 

K, 

1 — K 

3' 

3 

-1 

0 

1 — K 

K 

4 

4 

0 

1 

-1 

-1 

5 

5 

1 

-1 

0 

0 


Table 3: 


The character table of the group, where k 


_ 1+75 
2 


where R represents any irreducible representation of A^, and 4i, 42 , 5i and 62 stand for the 
two 4 and two 5 representations that appear in the Kronecker products. 

We now list the Clebsch-Gordan coefficients for our basis. We use the notation ctj (Pi) 
to denote the elements of the hrst (second) representation. The subscript ”S” (”A”) refers 
to symmetric (antisymmetric) combinations. 


3 0 3 = Is © 3 a © 5s 

3' © 3' = Is © 3)^ © 5s 

3 © 3' = 4 © 5 

Is : Q;i/3i + (+ 2/53 + ci3/^2 
/Cl2/^3 ~ Cl3/^2\ 

3a : 1 cii/^s — ci2/3i j 
\ci3/3i ~ 

/2ail3i — a2P3 — ci3/32\ 

— V^(cIi/ 32 + Cls/^i) 

5s ; v/6ci2/32 

VGas/Ss 

\ + as/Si) J 

Is : OiiPi + (^2^3 + CI 3/32 
/ <^2^3 — Cl3/^2\ 

3A • ( C1i/ 32 — Cl2/3l j 
\as/3i - aiPsJ 

/2ai(3i — 0:2/33 — 03/32^ 
VGasPs 

5s : —V^(cii/ 32 + ci2/3i) 

— V^(oi/33 + 03/3i) 

V ^ 602/^2 / 

/ V^Cl2/3l + C13/32 \ 

^ . —V^Oii/32 — OisPs 

— a/20i/ 33 — 02/32 
\ V^Cl3/^l + C12/33 / 

/ \/3oi/3i \ 

Os/^i — V^Os/^s 

5 : Oi/32 — 

OilPs ~ V^Cl2/32 
\03/3 i — \/2a2/33/ 


3©4 = 3'©4©5 

3'©4 = 3©4©5 

/ — 7 ^( 02/54 + Cl3/^l) \ 

3' ; 1 y/2ai/32 — 02/3i + | 

/ aijSi — \/2a3P2 \ 

^ _ —Oi/ 32 — V^Cl2/3l 

Cll/33 + V^Cl3/34 
\—O 1/34 + v^Os/^s/ 

/ ^ 6 ( 02/34 - a3/5i) \ 

2^/201/31 + 203/^2 

5 ; —V^oi/ 32 + Os/^i + 3 o3/33 

'\/20i/33 — 302/32 — 03/34 
\ —2v^oi/34 — 202/^3 / 

/ — 7 ^( 02/53 + 03 /^ 2 ) \ 
3 ; j v^cii/3i + O 2/34 — O 3/33 j 
\V^CIi/54 ~ + 03/3i/ 

/ oi/3i + ^ 

^ . cii/32 — 

—aiPs + \/2a2Pi 
X—aiPi — v^os/^s/ 

/ 7 ^( 02/53 — 03 /^ 2 ) ^ 

V^Oi/^i — 302/34 — 03/^3 

5 ; 2v^Oi/32 + 203/^4 

—2^/201/33 — 2o2/3i 
7 — 7 ^ 01/34 + 02/^2 + Sa^Pi/ 
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305 = 3© 3'©4©5 


3'05 = 3©3'©4©5 


3' 


4 ; 


—20Lij3i + '\/3ci2/^5 + '\/3Q^3/^2\ 

\/^aif32 + Q;2/5i — 1 

a/Soi/Ss - A/6tt2/54 + tts/Si / 
+ 0:2/35 + 03/32 \ 
01/33 — \/ 2 a 2 l 32 — \/ 2 a‘il 3 i 1 
01/34 — v^q;2/33 — 
f 2 y/ 2 aiP 2 ~ VQo' 2 Pi + O3/33 ^ 
—'\/2oi/33 + 20 : 2/32 — 803/34 
v / 2 oi /34 + 802/33 — 203/35 
2v^Oi/35 — O 2/34 + y/Qa^j3iJ 


'\/3oi/3i + O 2/34 + O 3/33 \ 
oi/32 — v^ 02/35 — v/203/34 1 
O 1/35 — V^02/33 — v/ 203 / 32 / 

—20i/3i + V^02/34 + V^03/33\ 
3^ ; ( v^Oi /33 + 02/3i — y/Qa^jS^ 1 

a/30i/34 - a/602/32 + 03/3i / 
/ y/2aiP2 + 802/35 — 203/34 ^ 
2v^oi/33 — v/6o2/3i + 03/35 
-2^/201/34 — 02/32 + '\/6o3/3i 

V ~v^oi/ 35 +202/33 — 803/32 y 



V^(o2/35 — O3/32) ^ 



v/3(o2/34 — O3/33) ^ 



— Oi/32 — V^02/3i — V^03/33 



2oi/32 + y / 2 asl 34 


5 ; 

— 20i/33 — v^02/32 


5 ; 

—01/33 — v^02/3i — y / 2 asP ^ 



2oi/34 + 



O1/34 + '\/2o2/32 + '\/3o3/3i 



\ 01/35 + V^02/34 + y/SasPi / 



\ —20i/35 — V^02/33 / 



4 0 4 = Is © 3 a © 3^ © 4s © 5s 

405 = 3© 3'©4© 5i© 52 

Is : O1/34 + 02/33 + O3/32 + 04/3 i 

/—O1/34 + 02/33 — 03/32 + 04/3i\ 

3a • 1 V^(02/34 — O4/32) 1 

\ V^(oi/33 — 03/3i) / 

/ O1/34 + 02/33 — 03/32 — 04/3i\ 

3a : 1 V^(o 3/34 — O4/33) j 

\ V^(Oi/32 — 02/3i) / 

/ 02/34 + 03/33 + 04/32'^ 

^ _ oi/3i + 03/34 + 04/33 
^ Oi/32 + 02/3i + 04/34 

\Oi/33 + 02/32 + 03/3i y 

/ V^(Oi/ 34 — O2/33 — O3/32 + 04/3l)^ 
—V^02/34 + 2y/2asl3s — v^04/32 

5s : —2^/201/31 + '\/2 o3/34 + y/2a4l3^ 

\/2oi/32 © \/2o2/3i — 2v^04/34 
V^Oi/33 + 2v^02/32 — y/^a^PiJ 

/ 

4 ; 

V 

51 ; 

52 : 

2y/2aij3^ — ^/2o2/34 + y/^a^Ps — 2y/2a4/32\ 

— V^Oi/3i + 2 o2/35 + 803/34 — Oi4j3‘4 1 

^ oi/34 — 802/33 — 2o3/32 + '\/6o4/3i y 

'^V^Oi/35 + 2^/202/34 — 2^/203/33 — V^04/32\ 
8oi/32 — '\/6o2/3i — 03/35 + 204/34 1 

y —20 i/33 + 02/32 + '\/603/3i — 804/35 y 

v/30i/3i — V^02/35 + V^03/34 — 2^/204/33 \ 

— V^Oi/32 — V^02/3i + 2^/203/35 + V^04/34 
V^Oi/33 + 2^/202/32 — V^03/3i — ^/204/35 

— 2V^Oi/ 34 + V^02/33 — V^03/32 + V^04/3iy 
(yf2a\fi^ — v^02/34 — y/2a^l3z + v^04/32\ 

—V^oi/3i — v/303/34 — v/304/33 

V^Oi/32 + V^02/3i + V^03/35 
y/3a2l32 + ©^03/34 + y/3a4/3^ 

V —V^Oi/34 — V^02/33 — V^04/3i y 

( 2oi/35 + 4o2/34 + 4o3/33 + 204/32 \ 

4oi/3i + 2^/602/35 

—V^oi/32 + 2o2/3i — v^03/35 + 2^/604/34 
2'\/6oi/33 — v^02/32 + 2o3/3i — '\/6o4/35 

\ 2^/603/32 + 4o4/3i y 
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5 = 1, 


3'. 


is,i 


4a,2 


>S,1 


>S,2 


4s, 


3' 


1 • 


Is : Oi/^i + 02/35 + 03/34 + O4/33 + O5/32 

( O2/35 + 2O3/34 — 2O4/33 — 05/32 
— \/3oi/32 + \/3o2/3i + \/2o3/35 — \/2a^(5^ 

'v/3oi/35 + a/2o2/34 — a/204/32 — 'v/305/3i 
/ 202/35 — 03/34 + 04/33 — 205/32 \ 

: I \/ 3 oi /33 — \/ 3 o 3 / 3 i + \/ 2 o 4/35 — \/2 a ^[3 i j 
\/3oi/34 + \/2o2/33 — \/2o 3/32 + \/3o4/3iy 
/^3a/2oi/32 + 3a/2o2/3i — \/3o3/35 + 4-\/3o4/34 — •v^05/33^ 
31 ^/ 201/33 + 4\/3o2/32 + 3\/2o3/3i — \/3o4/35 — \/3o5/34 
3a/2oi/34 — 1/302/33 — 1/303/32 + 3i/2o4/3i + 4i/3o5/35 
\ 3 i / 2 oi /35 — 1/302/34 + 41/303/33 — 1/304/32 + 3 i / 2 o 5 / 3 i / 
/ i/2oi/32 — \/2o2/3i + \/3o3/35 — \ 

—1/201/33 + 1/203/31 + 1/304/35 — 1/305/34 
— \/2oi/34 — \/3o2/33 + \/3o3/32 + \/2o4/3i 
\ 1/2O1/35 — \/302/34 + \/304/32 — \/205/3i / 
/2oi/3i + 02/35 — 203/34 — 204/33 + 05 / 32 ^ 

Oi/32 + 02/3i + 1/603/35 + 1/605/33 
—2oi/33 + \/6o2/32 — 2o3/3i 
—2oi/34 — 2o4/3i + y / Qoi ^/ 3 ^ 

O1/35 + 1/6O2/34 + 1/6O4/32 + 05 / 3 i 
/2oi/3i — 202/35 + O 3/34 + O 4/33 — 205 / 32 ^ 

—2oi/32 — 2o2/3i 4“ \/6o4/34 
O1/33 + 03/3i + 1/604/35 + 1/605/34 
O1/34 + 1/6O2/33 + 1/6O3/32 + 04/3i 
—2oi/35 + — 2o5/3i 


4a, 


5s,: 


V 


/ 


>S,2 


V 


/ 


B Lepton flavor mixing from As family symmetry with¬ 
out CP 

In this section, we investigate the possible lepton mixing patterns which can be derived 
from only ^5 family symmetry without CP symmetry imposed. As usual, the three gener¬ 
ations of left-handed leptons are assigned to the triplet representation 3, and A 5 is broken 
into two different abelian subgroups Gi and Gi, in the charged lepton and neutrino sector 
respectively. The residual flavor symmetry G^ can only be a Z 2 or subgroup of A 5 since 
we assume neutrinos are Majorana particles here. In this approach, the PMNS matrix can 
be obtained by simply diagonalizing the representation matrices of the generators of Gi and 
Gu without resorting to the mass matrix [^|^|^. For G^ = K 4 and Gi is capable of 
distinguishing the three generations of charged lepton, i.e., the eigenvalues of the generators 
of Gi aren’t degenerate, the PMNS matrix would be completely hxed up to row and column 
permutations. However, only one column would be hxed by the remnant havor symmetries 
Gi and Gy in case of G^ = Z 2 . In the following, the scenario of Gi = Z 2 and G^ = K 4 shall 
be discussed as well, and one row would be hxed instead. It is noteworthy that two pairs of 
subgroups {Gi,G^) and (G/G'j,) lead to the same result for the PMNS matrix, if they are 
conjugate under an element of the A 5 group. 
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B.l 


G. = Ka 


From Appendix]^ we know that Gi can be a Z3, Z5 or 7^4 subgronp of A5. In case of 
G; = Z5, all 6 X 5 = 30 possible combinations of Gi and G^, give rise to the same mixing 
matrix 


UpMNS— 


J 1 0 \ 


2F5k 


\ Y 2^/5K 


1 

\/5k 

K 

2 V 5 

K 

2\/5 


1 

F2 

1 

F2 


= Ugr 1 


(B.l) 


which is the well-known golden ratio mixing pattern. The mixing angles are determined to be 
sin^ 012 = (3 — k) /5 ~ 0 . 276 , sin^ 023 = 1/2 and sin^ 0i3 = 0 . Obviously 0i3 should acquire 
moderate corrections to accommodate the measured non-vanishing value of the reactor angle 
although 012 and 023 are in the experimentally favored 3 a ranges [^. 

In case of Gi = Z3, we hnd two mixing patterns can be obtained. For the representative 
symmetries Gi = and G^, = s,tst elements of Gi and G^ generate 

an A4 subgroup instead of the full flavor symmetry group A5. The resulting mixing matrix 
is given by the familiar democratic mixing in which all elements have the same absolute 
value [^, i.e., 

1 1 1 

TTl 

1 —6 3 

1 


UpMNS — 


27ri 

6 3 


= U, 


DM 


—e 3 


27ri 

e 3 


(B. 2 ) 


The mixing angles are sin^ 0i2 = sin^ 023 = 1/2 and sin^0i3 = 1 / 3 . Large corrections to 0i2 
and 013 are needed to be compatible with experimental data. For another representative 
symmetries Gi = and G^ = st st ) ^ parent group A5 can be generated by 

Gi and Gj,. The mixing matrix is found to be of the form: 


UpMNS — 

V6 



V2(l 


K, 


K 


K 



(B. 3 ) 


which leads to the following mixing angles: sin^ 0i2 = (2 - k) /3 ~ 0 . 127 , sin^ 023 = 1/2 and 


sin 013 = 0 . Notice that both 0i2 and 0i3 are outside of the 3 a ranges 5 . The same results 


have been obtained in Refs. 49,51 . For the last case of Gi = K4, where G^ and Gi are not 


the same Klein group, only one mixing pattern can be derived. 


K 

- 1 

-1 

1 — K 

K — 1 

K 

to 

{ 3 -k 



= Urc ■ 


(B. 4 ) 


and sin^ 013 = (2 — k) /4 ~ 0 . 0955 . Both 0i3 and 023 are too large to be acceptable. This 


mixing pattern has also been found in Ref. 49 . In summary, no mixing matrix in agreement 


with experimental data can be obtained if the full Klein symmetry is preserved by the 
neutrino mass matrix. In the following, we consider the situation with a single residual Z2 
flavor symmetry in the neutrino sector or in the charged lepton sector. 


B .2 Gi^ — Z2 or Gi — Z2 

In this case, only one column or one row of the PMNS matrix would be determined 
up to permutations and phases of its elements by the remnant flavor symmetries G/ and 
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Gi 

Gy 

Fixed column or row 


7 '^ 

yS 

I' nr 1 1 \T 

^ \l Gb' y/2VbG y/2VbK 

/ 

T'iST'^ST^ 

^2 


/ 

T^ST'^ST^S 

^2 


X 

T^st^S 

yS 


X 

yT^ST'^ST^ 

^2 

/ 1 — K K Hi \T 

'' \/3 ’ \/6’ 

X 

yT'i ST'^ ST^ S 
^2 

( K K— 1 Hi—1\T 

^ Vs’ V6 ’ Vb ^ 

X 

yT^ST^S 

yST'^ST^S 

^2 

( 1 1 1 \T 

''Vs’ Vs’ Vs’^ 

/ 

r^(ST^ST^S, TST*) 

yS 

^2 

f K, 1 K—1\T 

V 2 ’ 2 ’ 2 7 

/ 

j^iST^ST^S, TST'^) 

yTST^ 

^2 

(1,0,0)^ 

X 

z! 

j.iST'^ST^S, TST'^) 
-^4 

H Hi 1 Hi—1 \ 

1 2 ’ 2 ’ 2 7 

X 

yTST* 

^2 

r^iST^ST^S, TST'^) 
-^4 

(1,0,0) 

X 


Table 4: The possible form of one column (row) of the PMNS matrix determined by the residual flavor 
symmetry Gi, = Z 2 {Gi = Z 2 ) within the framework of flavor symmetry. The notation denotes 
that the relevant lepton mixing is compatible with the experimental data at 3a level |^. The notation “/f” 
implies the resulting mixing is not viable. 


Gy 34,52 . This method generally allows us to obtain relations between mixing parameters 
and a non-zero ^ 13 . We have scanned all independent combinations of Gi and Gy, and the 
corresponding explicit forms of the hxed column or row vector are presented in Table 
Comparing with the present 3a conhdence level ranges of the moduli of the elements of the 
leptonic mixing matrix 


\UpMNs\ 3 a — 


/ 0.789 ^ 0.853 
0.194 ^ 0.558 
\0.194 ^ 0.558 


0.501 ^ 0.594 
0.408 ^ 0.735 
0.408 ^ 0.735 


0.133 ^ 0.172 \ 
0.602 ^ 0.784 
0.602 ^ 0.784 / 


(B.5) 


we hnd that neither of the two possible row vectors can be accommodated by the data, and 
only four cases are viable. The remnant symmetries can be chosen to be (G;, Gy) = [Z^, Z^), 
T^ST^ST^\ ^ ^ST'^ST^s^ S',r 5 r witliout loss of generality, 


and the hxed column are 




\/ 2\/5 K' \/ 2\/bi- 


^{k, —1, K — 1)^ respectively. These column vectors can £t the hrst or the second column 
of the PMNS matrix. The resulting lepton mixing matrix can be obtained from Ugr, Urm 
and Urc by multiplying a unitary matrix U 23 or t /13 from the right-hand side with 


j 

/ cos 9 

0 

sin 9e \ 


(1 

0 

0 \ 

Un = 

0 

1 

0 

CO 

0 

cos6' 

sin6*e I 

1 

^ — sin 

0 

cos 9 ) 



- sin 9e^^ 

cos 6* / 


(B. 6 ) 


where 9 and 6 are real, and a arbitrary phase matrix in the right-hand side of 17i3 and U 23 
is omitted, since they can be absorbed into the Majorana phases which are not constrained 
by havor symmetry. The multiplication of U 13 {U 23 ) corresponds to performing a unitary 
linear transformation of the 1st (2nd) and 3rd columns. In the following, we shall discuss 
the predictions for the PMNS matrix and lepton mixing parameters in each case. 
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B.2.1 Gi = Zj,G, = Zi 


The lepton mixing matrix Upmns is predicted to have one column —j=^, —j==Y 

which coincides with the first column of the GR mixing. The other two columns should be 
orthogonal to it, and they can be obtained by making a unitary rotation of the 2nd and 3rd 
columns of Ugr- 


V a/5 


u. 


PMNS 


— UgrU23 — 


COS 9 

y/OK 


sin 9e 

VOK 




\f2^K 

1 




K a sine ^i(5 / K „• n^—iS i cos 9 

^cos9-^e ^^sm9e +^j 


(B.7) 


\ a/ 2 V 5 k 

This form of modification to the GR mixing has been discussed in a phenomenological way 


in Ref. 53 55 . Here we show that this mixing pattern can be naturally reproduced from 


the H 5 flavor symmetry. The mixing angles can be straightforwardly extracted as follows, 

3 — K, . n „ .n„ 2 cos^ 9 


sin 9 i 3 = 


sin^ 923 = - - 


sin 9 , sin 9 i 2 = 


3 + 2k + cos 26* ’ 


1 

2 


■\/3 + 4k sin 29 cos 5 


3 + 2k + cos 26* 

Wee see that the solar and reactor mixing angles are related by 

5 cos^ 6*12 cos^ 6*13 = 2 + K . 


(B.8) 


(B.9) 


Given the 3cr ranges 1.76 x 10 ^ < sin^ 6*13 < 2.95 x 10 ^ and 0.259 < sin^ 6 *i 2 < 0.359 from 
global analysis 1 ^) ^13 and 6*12 are further constrained to be in the intervals of 1.76 x 10 ^ < 
sin^ 6*13 < 2.35 x 10“^ and 0.259 < sin^ 6*12 < 0.263 by this correlation. The well-known 


Jarlskog invariant Jcp 56 , which measures the size of the GP violation, is written as 


Jgp = — ^ - sin 26* sin 6 . 

20 

The Dirac GP violating phase 6 cp is expressed in terms of 6 * and 6 as 


sin (5cp = 


1 / 2(3 -|- 2k -|- cos 26*)sign(sin 26*) sin 6 


\/ 4(3 -|- 2k) cos 26* -|- (7 -|- 8k) (3 -|- cos 46*) — 4(3 -|- 4k) cos 26 sin^ 26* 


(B.IO) 


(B.ll) 


In order to see how well the lepton mixing angles can be described by this mixing pattern 
and its prediction for 6cp, we perform a numerical analysis. The free parameters 6 * and 6 
are scattered in their whole allowed ranges of 0 < 6 * < 27r and 0 < 5 < 27r. The correlations 
and the possible values of the mixing parameters are plotted in Fig. Furthermore, the 
experimental data of three mixing angles 6 * 12 , 6*13 and 6*23 at 3cr level are considered, 
accordingly the allowed values of the mixing parameters would generically be constrained 
in small regions. Here and herafter, we perform numerical analysis and present results only 
for normal ordering neutrino mass spectrum. The results would change a little bit for the 
inverted ordering case. From Fig. we can read that sin^ 6*12 is predicted to be around 0.26, 
any value of 6*23 within the 3a range can be achieved and 6cp is restricted in the range of 
[0.990,2.152] U [4.131,5.293]. Recalling that if both A 3 family symmetry and generalized 


GP are imposed, as discussed in section 3.1, the parameter 6 can only be 7 r /2 (case H) 
rather than free. Note that case I is not viable. As a consequence, the Dirac GP 6 cp would 
be maximal. Therefore we conclude that the generalized GP symmetry is a quite effective 
method of predicting the GP violating phases. 
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Figure 6 : Predictions for the mixing parameters sin^ 0 i 2 , sin^ 6 * 23 , Jcp and 6cp with respect to sin 6*13 when 
the remnant flavor symmetries are Gi = and Gi, = Zf. The corresponding PMNS matrix is given by 
Eq. (B.7). The pink regions denote the possible values of the parameters when both 0 and 5 freely vary in 
the whole region of [0, 27r]. The dark green areas represent the regions allowed by the current experimental 
data for three neutrino mixing angles at 3a level [^. The red pentagrams refer to the best fitting values of 
case II discussed in section [XT] after the generalized CP is imposed. 


B.2.2 Gi = Zj,G, 


^T^ST'^ST^ 


In this case, one column of Upmns is determined to be (y 

is exactly the second column of the GR mixing. The corresponding 
obtained from Ugr by multiplying f/13 from right-hand side, 



PMNS matrix can be 


UpMNS — UgrUis — 


-^cos9 

cos 6 I sin 6 iS 

cos 9 _ sin 9 ^iS 

\ \/2vTk ^ 


The lepton mixing parameters read 

sin^ 6*13 = ^ sin^ 9 , sin^ 6*12 = 


1 

\/5k 

K 

2 V 5 

K 

2V5 


- ,/^sin0 e-*'^ \ 

Y v5 

_ cos 9 I _sin_0_ — 2 ^ 
cos 9 I sin^ 


3 -b /t -b (1 -b /t) cos 29 ’ 


(B.12) 
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Figure 7: Predictions for the mixing parameters sim^i2, sin^^23, Jcp 5cp with respect to sin ^13 
when the remnant flavor symmetries are Gi = and Gi, = . The corresponding PMNS matrix 

is given by Eq. (B.12|. The pink regions denote the possible values of the parameters when both 9 and 6 
freely vary in the whole region of [0, 27r]. The dark green areas represent the regions allowed by the current 
experimental data for three neutrino mixing angles at 3a level [^. The red and yellow pentagrams denote 
the best fitting values of case III and case IV discussed in section where the generalized CP symmetry 
is considered. Notice that the red pentagrams almost coincides with the yellow one in the first panel, since 
the best fitting values of sin^ 0 i 2 and sindia are nearly the same in case III and case IV. 


sin^ 023 = 
sin 6c p = 


\/2 + K sin 29 cos 6 


a/2 + K ■ OQ ■ X 

^CP = —^—sin 20 sin 0, 


2 3 +/t + (1 + k) cos20’ 20 

^2{2 + k){?>k-2 + K COS 20) sign(sin 20) sin 6 
a/(13 + An) (3 + cos 40) + 4(7 + Ok) cos 20 — 20 sin^ 20 cos 26 


We have a relation between 0i2 and 0i3, 


(B,13) 


5 sin^ 012 cos^ 013 = 3 —/t. (B.14) 

The solar mixing angle 0i2 is restricted by the observed value of 0i3 such as 0.281 < sin^ 0i2 < 
0.285 which is in the 3a range [^. We display the allowed regions of the mixing angles, Jcp 
and 6 cp in Fig. No dependence of 6 cp on sin0i3 is observed, and 6 cp can take any value 
in the whole range of [0,27r]. However, 6 cp can only be conserved or maximally broken 
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if generalized CP is considered, as shown in section 3.2 Note that the mixing pattern in 


Eq. (B.12) has been discussed in Ref. 53,55 


B.2.3 Gi = Z] 


\ G. = Zi 


1 T 

The chosen remnant symmetry leads to a trimaximal column (1,1,1) , and Upmns 
takes the form 


cos 9 


AS 


sin 6* 


t/, 


PMNS 


= UncUi-i = —= \ cos6'+ 6*^3+'^) sin6' 


73 


1 cos 9 + e sin 6* 

1 sin^ — cos6' 


—e’s’cos^ — sin^ 1 e^3* cos^ — e**^3 '^)sin6' 


(B.15) 

Such a mixing pattern as a minimal modihcation to the tri-bimaximal has been widely 
discussed in the literature 53,55 IHZl. and it can also be naturally reproduced from simple 

The predictions for the lepton mixing 


flavor symmetries 23,58 and 5*4 |20[|24 
parameters are given by 

sin^ 6*13 = x(l + sin 26* cos 5), 


sin^ 6*10 = 


sin^ 6*23 = V 


sin 6 c p = 


■y/S sin 26* sin 5 
4 — 2 sin 29 cos 6 ’ 


Jcp — 


2 — sin 29 cos 6 ’ 
cos 26* 


673 ’ 


-72 cos 26*(2 — sin 26* cos <5) 


da — sin 26* cos 5) (5 + 3 cos 46* + 2 sin^ 26* cos 36) 
As expected, the following relation is fulfilled, 

3 sin^ 6*12 cos^ 6*13 = 1, 


(B.16) 


(B.17) 


which generically holds true for trimaximal mixing. Inserting the experimental bound of 
^13 we obtain 0.339 < sin^ 6*12 < 0.343. A numerical analysis similar to previous cases is 
performed, as shown in Fig. We see that no prediction for 6 cp can be made. Recalling 
that 6 cp would be constrained to be maximal by generalized CP symmetry discussed in 
section 13.31 


B.2.4 


Ui — A 4 , (jr, 


Z 


s 

2 


One column is hxed to be ^ —1^ k, — 1)^ in this case, and it can only be the hrst 

column of the PMNS matrix in order to be consistent with the experimental data. As a 
result, UpMNS is of the form 


UpMNS — UrcU23 


( K — COS 9 — sin 9e^^ 
— 1 — COS0 — «:sin6*e*'^ 

K — 1 K cos 6* — sin 6*6**^ 


cos 9 — sin 9e~^^ 
K cos 9 — sin 9e~^^ 
cos 9 + K sin 6*6“**^ 


(B.18) 


Then the three mixing angles read 


sin^ 6*13 = 



(75 — cos 26* — 2 sin 26* cos 5), 


sin^ 6*12 = 
sin^ 6*23 = 


3 — K + (k — 1) (cos 26* + 2 sin 26* cos 5) 

5 + K + {k — 1) (cos 26* + 2 sin 26* cos 5) ’ 
3 + 75 cos 26* — 2 sin 26* cos 6 
5 + K + (k — 1) (cos 26* + 2 sin 26* cos 5) 


(B.19) 
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Figure 8: Predictions for the mixing parameters sim ^ 12 , sin^ Q- 
the remnant flavor symmetries are Gi = ^ and G 


23 , Jcp and (5cp with respect to sin 0 i 3 when 
The corresponding PMNS matrix 


_ ^ST^ST-'S 


is given by Eq. (B.15I. The pink regions denote the possible values of the parameters when both 9 and S 
freely vary in the whole region of [0, 27r]. The dark green areas represent the regions allowed by the current 
experimental data for three neutrino mixing angles at Scr level [^. The red pentagrams refer to the best 
fitting values of case V discussed in section [331 after the generalized CP is imposed. 


A relation between 612 and 6*13 follows immediately 

4 cos^ 6*12 cos^ 013 = 1 + K. (B.20) 

The solar mixing angle is predicted as 0.326 < sin^ 0i2 < 0.334 which is in the experimental 
3a bound [^. The Jarlskog invariant Jcp is given by 

Jcp = ——sin 20 sin 5 . (B.21) 

16 

The Dirac CP violating phase 5cp is 

—y/2K — 3(6k + 1 + cos 20 + 2 sin 20 cos 6 ) sin 20 sin 6 
sin dcp = , = 

Y [5 — (cos 20 + 2 sin 20 cos <5)^] (\/5 — cos 20 + 2 sin 20 cos 5) (3 + \/5 cos 20 — 2 sin 20 cos 5) 

(B.22) 
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sin6'i3 



Figure 9: Predictions for the mixing parameters sin^ 0i2, sin^ 6*23, Jcp and 6cp with respect to sin 6*13 when 
the remnant flavor symmetries are G/ = S,tst ) corresponding PMNS matrix 

is given by Eq. (B.I 8 I. The pink regions denote the possible values of the parameters when both 9 and 6 
freely vary in the whole region of [0,27r]. The dark green areas represent the regions allowed by the current 
experimental data for three neutrino mixing angles at 3 ct level [^. The red pentagrams refer to the best 


htting values of case VII with 023 ( 06 /) < 45° discussed in section 3.4 after the generalized CP is imposed. 


The numerical results are displayed in Fig. We see that 6 cp is predicted to be in the range 
of [0,1.043] U [5.240, 27r], and the atmospheric mixing angle 6*23 mostly is less than 45° (i.e., 
in the hrst octant) in order to be compatible with experimental data of 6^13. The scenario of 
623 in the second octant can be achieved, if the second and third rows of the PMNS matrix 
in Eq. (B.18) are exchanged. Then the predictions for the solar and reactor mixing angles 
in Eq. (B.19) remain, 6 cp becomes tt + 6 cp, and 623 becomes 7 i /2 — 6*23. Consequently both 
Jcp and sin(5cp change into their opposite, and the expression of sin^ 6^23 in Eq. (B.19) is 


replaced by 


sm 


^23 — 


k(\/ 5 — cos 29 + 2 sin 29 cos 5) 

5 + K + (k — 1) (cos 29 + 2 sin 29 cos 5) 


(B.23) 


The predictions for sin^6'23 and 6 cp versus sin^is are shown in Fig. 10 As expected, ^23 is 


really larger than 45° to accommodate the measured values of 6*13, and the CP phase 6 cp is 
in the range of [2.099,4.185]. Notice that generalized CP would constrain 6 cp to be trivial. 
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Figure 10: The correlations of sin^ 023 and 5cp with respect to sin0i3, where the PMNS matrix arises 
from an exchange of the second and third rows in the pattern in Eq.(B.18|. The pink regions denote the 
possible values of the parameters when both 0 and 5 freely vary in the whole region of [0,27r]. The dark 
green areas represent the regions allowed by the current experimental data for three neutrino mixing angles 
at 3(7 leve l [^ . The red pentagrams refer to the best fitting values of case VII with 023 (^&/) > 45° discussed 
in section [3^ after the generalized CP is imposed. 


as studied in section 3.4 In summary, if a single Z 2 subgroup of the flavor symmetry 


is preserved by the neutrino mass matrix, only one column of the PMNS matrix can be 
determined and agreement with experimental data can be achieved. However, the Majorana 
phases cannot be predicted by flavor symmetry, and the Dirac phase 5cp is constrained very 
weakly. On the other hand, if we extend the family symmetry to include the generalized 
CP, 5cp is predicted to be trivial or maximal and Majorana phases are trivial. 
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